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We study the replica field theory which describes the pinning of elastic manifolds of arbitrary internal dimen- 
sion d in a random potential, with the aim of bridging the gap between mean field and renormalization theory. 
The full effective action is computed exactly in the limit of large embedding space dimension N. The second 
cumulant of the renormalized disorder obeys a closed self-consistent equation. It is used to derive a Functional 
Renormalization Group (FRG) equation valid in any dimension d, which correctly matches the Balents Fisher 
result to first order in e — 4 — d. We analyze in detail the solutions of the large-TV FRG for both long-range and 
short-range disorder, at zero and finite temperature. We find consistent agreement with the results of Mezard 
Parisi (MP) from the Gaussian variational method (GVM) in the case where full replica symmetry breaking 
(RSB) holds there. We prove that the cusplike non-analyticity in the large TV FRG appears at a finite scale, 
corresponding to the instability of the replica symmetric solution of MP. We show that the FRG exactly repro- 
duces, for any disorder correlator and with no need to invoke Parisi's spontaneous RSB, the non-trivial result 
of the GVM for small overlap. A formula is found yielding the complete RSB solution for all overlaps. Since 
our saddle-point equations for the effective action contain both the MP equations and the FRG, it can be used to 
describe the crossover from FRG to RSB. A qualitative analysis of this crossover is given, as well as a compar- 
ison with previous attempts to relate FRG to GVM. Finally, we discuss applications to other problems and new 
perspectives. 



I. INTRODUCTION 

Elastic objects pinned by a quenched random potential are a 
relevant model for many experimental systems. It describes 
interfaces in magnets (jj, |2| which experience either short- 
range disorder (random bond), or long range (random field) 
disorder, the contact line of a liquid wetting a rough substrate 
O10I- vortex lines in superconductors A 0,0 HI- It also pro- 
vides powerful analogies, via mode coupling theory, to com- 
plex systems such as structural glasses [9]. One important 
observable is the roughness exponent C, of the pinned mani- 
fold. 

From the theoretical side, this problem still offers consid- 
erable challenges. It is the simplest example of a class of 
disordered systems, including random field ma gnet s, where 
the so called dimensional reduction d OS HHITl flU Q 
renders conventional perturbation theory trivial and useless at 
zero temperature. The elastic object is usually parameterized 
by a N component vector u(x) in the embedding space R , 
and x £ K d is the coordinate in the internal space. Apart from 
the case of the directed polymer (DP) in 1 + 1 dimensions 
(d = 1, N = D, where some exact results were obtained 
Hi E d E El, analytical results are scarce. One im- 
portant challenge is to understand the DP for any N, due to 
its exact relation to the Kardar-Parisi-Zhang growth equation 
whose upper critical dimension is at pres ent not known, and 
even its very existence is debated 12(il2lll2^l23ll . 

Two main analytical approaches have been devised so far. 
Each succeeds in evading dimensional reduction, providing 
an interesting physical picture, but comes with its limitations. 
The first one is the mean field theory, the replica gaussian vari- 
ational method (GVM) [24] in the statics and the off equilib- 
rium dynamical version H25L 12611 . The GVM approximates the 



replica measure by a replica symmetry broken (RSB) gaus- 
sian, equivalently, the Gibbs measure for u as a random su- 
perposition of gaussians [24], and is argued to be exact for 
/V = oo. It yields Flory values for the exponent £. As for spin 
glasses, computing the next order corrections (i.e. in 1/JV) at 
the RSB saddle point is very arduous 11271 l28l 12911 . One may 
question whether it is the most promising route, since it is as 
yet unclear whether the huge degeneracy of states encoded in 
the Parisi RSB is relevant to describe finite N. There seems 
to be some agreement that this type of RSB does not occur 
for low d and N. Certainly, in the simpler but still-non trivial 
d = limit, Parisi type RSB found in the GVM should exist 
only at N — oo, apart from the interesting so-called marginal 
case of logarithmic correlations H30I1 . For the DP, another ex- 
actly solvable mean field limit is the Cayley tree and there too 
it is not clear how to meaningfully expand around that limit 

immin. 

The second main analytical method is the functional renor- 
malization group (FRG) which performs a dimensional ex- 
pansion around d = 4 and was originally developped only to 
one loop, within a Wilson scheme l a l34Tl35ll36ll . Its aim is to 
include fluctuations, neglected in the mean-field approaches. 
There too, the dynamics has been investi- 

gated. The FRG follows the second cumulant of the random 
potential R(u) under coarse graining, a full function since the 
field is dimensionless in d = 4. It was found that R(u) be- 
comes non-analytic already in the 1-loop equation at T = 
after a finite renormalization, at the Larkin scale. 

Both methods circumvent dimensional reduction by pro- 
viding a mechanism which is non-perturbative in the bare 
disorder. The GVM evades DR thanks to the RSB saddle 
point. The FRG escapes via the generation of a cusp-like 
non-analyticity in R"(u) at u = 0. Indeed, while the bare 
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disorder correlator is an analytic function, FRG fixed points 
for the renormalized R(u), perturbative in e = 4 — d, are 
found only in the space of non-analytic functions, and subject 
to the condition that the resulting exponent ( is non-trivial. 
Both methods are disconcertingly different in spirit and it is 
an outstanding question in the theory of disordered systems 
how to compare and reconcile them. Comparisons were made 
between some predictions of the 1 -loop FRG and of the GVM 
161 13611 . Balents and Fisher obtained the 1-loop FRG equation 
for any TV restricted to 0(e), and found that its solution repro- 
duces the Flory value of £ for LR disorder, but yields subtle 
corrections for SR disorder, exponential in TV. 

Physically both methods capture the metastable states be- 
yond the Larkin scale L r and it is tempting to compare how 
they describe them. In 1 42] a coarse grained random poten- 
tial was defined and it was found within the GVM that its 
correlator mimics the one in the FRG, exhibiting some non- 
analyticity which was interpreted in terms of shock-like sin- 
gularities in the coarse grained disorder. Unfortunately, this 
analogy was demonstrated only around the Larkin scale, while 
a quantitative and more general connection able to reach per- 
turbatively the true large scale behaviour, as is achieved in the 
field theoretic FRG, is still missing. 

The need for a study of the FRG at large TV is all the 
more pressing since we have developped systematic higher 
loop approaches within the e-expansion j43l l44l l45il . Within 
these studies, we have found that higher loop FRG equations 
for R(u) at ti ^ contain non-trivial, potentially ambigu- 
ous "anomalous terms" involving the non-analytic structure 
of R(u) at u = 0. We have proposed a solution to lift these 
ambiguities in the statics at two loops i43ll44l l45ll . Since the 
large-TV limit allows in principle to handle higher-loop cor- 
rections (i.e. to treat any e) it should be useful to understand 
the many-loop structure of the field theory. Stated differently, 
we want to understand which physical quantity precisely does 
the FRG computes? Finally, developping a systematic 1/TV 
expansion within the FRG for any d should provide a novel 
handle to attack problems such as KPZ, maybe avoiding the 
need for spontaneous RSB altogether if it proves to be non- 
essential. 

The aim of this paper is to study the FRG at large TV. For 
this purpose we first perform an exact calculation of the effec- 
tive action T[u] of the replicated field theory at large TV. Its 
value for a uniform mode and further expansion in cumulants 
yields a definition of the renormalized disorder consistent with 
field theoretic approaches. The second disorder cumulant is 
found to obey a closed self-consistent equation. All higher cu- 
mulants can be constructed recursively from the lower ones. 
It can be easily inverted below the Larkin scale and there the 
solution is analytic and corresponds to the replica symmetric 
solution of MR Varying with respect to an infrared scale, here 
the mass, we obtain the FRG /3-function in any d at dominant 
order, TV = oo. The continuation beyond the Larkin scale 
is remarkably easier to perform on the resulting FRG equa- 
tion. Its solution reveals that the FRG exactly reproduces the 
non-trivial result of the GVM with full RSB for small over- 
lap. We also give a formula which yields the complete RSB 
solution for all overlaps. At no point in our derivation Parisi- 



RSB is invoked, as replica symmetry is broken explicitly here. 
Since our saddle point equations for the effective action con- 
tain both the MP equations and the FRG, it can be used to 
describe the crossover from FRG to RSB. A qualitative anal- 
ysis of this crossover is given, as well as a comparison with 
previous attempts to relate the FRG to the GVM 1 42]. Finally, 
applications to other problems and new perspectives are dis- 
cussed. A short version of this work has appeared in |46J]. In 
a related paper 14711 . we give all details of the calculation of 
the 0(1/ N) corrections, with the aim of understanding finite 
but large TV. 

The outline of the paper is as follows. In SectionllTlwe de- 
fine the model, the effective action and its physical interpre- 
tation. In Section|n]we compute the effective action at large 
TV, using the saddle point method and perform a cumulant 
expansion (Section IFvl . A graphical interpretation is given 
in Section|V] In Section lvTl we establish the FRG equation at 
large TV (the /3-function of the theory). Then in Section lVTTI we 
perfom a detailed analysis of the FRG equation for a specific 
class of disorder correlators, both below and above the Larkin 
scale. In Section lVllll we compare the FRG with the MP so- 
lution using RSB. First we recall the MP approach and find 
agreement with the predictions of the FRG calculation. Next 
we extend these results to an arbitrary disorder correlator for 
which the GVM gives full RSB. Finally we discuss the phys- 
ical interpretation and compare our approach with the one of 
Ref. 1 42]. Section lDTI presents the conclusion. The appendices 
contain several generalizations, the calculation of the third and 
fourth disorder cumulant, finite temperature fixed points, and 
an analysis and comparison with the effective action in more 
conventional field theories. 

II. MODEL AND PROGRAM 

A. Model and large-TV limit 

We consider the general model for an elastic manifold of in- 
ternal dimension d embedded in a space of dimension TV. 
The position of the manifold in the embedding space is de- 
scribed by a single valued displacement field u(x), where x 
belongs to the internal space and u is a TV component vector 
which belongs to the embedding space. (Its components it 1 , 
i = 1 , . . . , TV, are specified below only when strictly neces- 
sary.) A well studied example is that of an interface (e.g. a do- 
main wall in a magnet) where d = 2 and TV = 1. There u(x) 
denotes the height of the interface. Other examples are the di- 
rected polymer {d = 1) in a TV dimensional space, which can 
be mapped to the TV-dimensional Burgers and Kardar-Parisi- 
Zhang (KPZ) equations 14811 . or a vortex lattice in the absence 
of dislocations described by d = 3, TV = 2, where u(x) is 
there the deformation from the ideal crystal 0,0- 

We will study here the equilibrium statistical mechanics of 
such an elastic manifold in presence of quenched disorder, 
modeled by a random potential V(x, u(x)). It is described, 
in a given realization of the random potential, by the partition 
function 




-Hv[u]/T 
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where 

H v [u] 



-C(q)- L u(-q)-u(q) + T / V(x,u(x)) 



(2.2) 

consists out of an elastic energy (expressed here in Fourier 
space and taken to be isotropic), and of a pinning energy due 
to disorder. Here and below we denote 



d d q 



d d x 



(2.3) 



and u ■ v = Y^iLi u%v% ■ Throughout, square brackets as e.g. 
in A[u] denote a functional, here A of the field u a (x), while 
parenthesis as in A{u) denote functions. 

A convenient form for the inverse bare propagator, used be- 
low, is: 



T 



(2.4) 



where T is the temperature and the elastic constant is set to 
unity by a choice of units. The role of the additional mass 
term m will be discussed below. An additional small scale 
(ultraviolet, UV) cutoff A is implied here and will be made 
explicit when needed. 

This model is highly non-trivial and, apart from the cases 
of N = 1 and d = 0,1 , very few exact results are known 
E E E3 E d Eai- To obtain exact results for large 
embedding space N — > oo, we need to consider a fully 
isotropic version of the model with 0(N) symmetry such that 
the model remains non-trivial in that limit. As in standard 
large- N treatment (as for instance of the <fi 4 O(N) model) 
one defines the rescaled field 



v(x) 



u(x) 



N 



(2.5) 



We will freely switch from one to the other in the following. 
One also chooses the distribution of the random potential to 
be O(N) rotationally invariant. It can be parameterized by its 
set of connected cumulants, of the form 



V(x, u)V(x', it') = R(\u - u'\)S d {x - x 1 ) 

= NB({v-v') 2 )5 d (x-x') (2.6) 



V{xi,ui) . . . V(x p ,u p ) 

= NS xl ,... tXp (-iysM(v 1 ,...,v p ), p>3 (2.7) 



$xi,...,x p ■= Y\_d d (xi - Xi) 



(2.8) 



This adequately models the case of uncorrelated (or short- 
range correlated) disorder in the internal space, studied here. 
The second cumulant, which plays the central role, is thus de- 
fined in terms of a function B(z). The higher cumulants are 
not strictly necessary in the bare model, but they appear, as 
we will see, under coarse graining. The distribution of dis- 
order being translationally invariant, these functions satisfy 



S (p){ 



Vi 



+ v) = S^(vi, . . . , v p ) for any v. The 



model studied here is thus a slight generalization of the model 
studied by Mezard and Parisi 1 24], henceforth also referred to 
as MP, in the same limit. 

Although we will consider the general case, it is useful, as 
in MP 1 24] to define two sets of simple models for which 
more specific results will be given. These are, respectively, 
the gaussian, short-range (SR) disorder, correlator 



B{z) = ge~ z 
and the power-law correlations 

B[z) = ey A i) (a2 + z)1 ^ 



(2.9) 



(2.10) 



which, for infinite N always corresponds to long-range (LR) 
disorder, a different universality class, as we will see below. 
For finite N, the long-range disorder corresponds, at the bare 
level, to 7 < 1 + A/2; but this is modified at the renormalized 
level, and the true frontier LR-SR for finite A is non-trivial. 

B. Program 

Having defined the model, and before turning to calculations, 
let us first outline what we aim at. All the considerations in the 
present section are valid for any A, but, since in the next sec- 
tion we will consider the large A limit explicitly, we already 
make apparent the rescalings. 

The model defined above has already been studied in MP 
I24B . One of the aims of this study was to compute the rough- 
ness exponent of the manifold, defined from the 2-point func- 
tion as 



((u(x) - u(x')) 2 ) ~A\x-x'\ 



'|2C 



(2.11) 



Besides the roughness exponent £, the amplitude A is also 
of interest whenever it is universal, as it is the case e.g. for 
long range disorder. To this aim the model was replicated 
(u — > u a ), averaged over disorder and self-consistent saddle 
point equations where derived for the 2-point function 



Gab(q) = (v a (q)v b (-q)) 



(2.12) 



This can always be done in a large- A limit, and is then solved 
via a RSB ansatz. 

Our goal is in a sense broader. We want to understand the 
full structure of the field theory, i.e. all correlation functions 
and not only the 2-point one. We will thus instead study the 
generating function of correlations as well as the effective ac- 
tion functional which yields the renormalized vertices. This 
program, defined here, will be carried out in the following 
sections explicitly for large A. In this article we will restrict 
ourselves to dominant order, but the aim is to understand large 
but finite N, including calculating of 1/A corrections. This is 
deferred to 14711 . 

1 . Effective action and field theory 

All physical observables for any N can be obtained from the 
replicated action in presence of a source, i.e. an external force 
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J a (x) acting on each replica: 



Z[J] = J Ij2?[u«]( 



■E,HvK]/T+J,E, J o (*)-Ua(a0 



(2.13) 

where u a (x), a = 1, . . . , n are the replicated fields (each one 
being an N component vector u a (x)). Differentiating with re- 
spect to the replicated source J a (x) in the limit n — > yields 
all correlation functions. The finite-n information is also in- 
teresting. For instance from 



Z := Z[J = 0] = exp(-nJT y /T) 



(2.14) 



one can retrieve the sample to sample distribution of the free 
energy Ty = — T In Zy, as was done e.g. in a finite size sys- 
tem for d = 1 1171 1.5011 . Thus, unless specified we will keep n 
arbitrary. 

One can explicitly perform the disorder average in (12. 1 3i : 



Z[J] = J ]]v[u a ]e- NS ^ 



.i] 



(2.15) 



+ [U( X (x)) - j a (x) ■ v a (x)} ,(2.16) 



S[u,j} = - / CiqY^i-q) ■ v a (q) 



where here Xab(x) = v a (x) ■ Vb(x) and here and below sum- 
mations over repeated replica indices are implicit. We have 
rescaled the source in a manner complementary to the field: 



J a (x)=VNj a {x). (2.17) 
We have also introduced the bare interaction potential 

U(X) = l^^B{Xab)--^^S(Xab,Xbc,Xca) + --- , 

ab abc 

(2.18) 

which is a function of a n by n replica matrix Xab and has a 
cumulant expansion in terms of sums with higher numbers of 
replicas. Because of translational symmetry and O(N) invari- 
ance it depends only on the matrix 



Xab ■— Xaa + Xbb — Xab — Xba 



(2.19) 



and the form of each cumulant is restricted. For instance one 
has SW(vi,v 2 ,v 3 ) = S((vi - v 2 ) 2 , (v 2 - v 3 ) 2 , (v 3 - vi) 2 ) 
etc.. The matrix potential U(x) can thus be considered as 
a convenient way to parameterize the disorder (here the bare 
disorder). 

The physical object which contains the information about 
the field theory at large scale is the effective action. It is the 
generating function of the 1 -particle irreducible diagrams and 
in conventional field theories its formal expansion in powers 
of the field yields the renormalized vertices. All correlation 
functions are then obtained simply as tree diagrams from these 
renormalized vertices. In particular it is known that within a 
d = 4 — e expansion at zero temperature to at least 2-loop 
order the theory can be renormalized (i.e. rendered UV finite 



and yielding universal results) by considering counter-terms 
only to the second cumulant. The latter is a function R(u), 
and can be viewed as the set of all coupling constants which 
simultaneously become marginal in d = 4. To probe renor- 
malizability to any number of loops, we want to compute the 
effective action from first principles. 

The effective action functional is defined as a Legendre 
transform: 



I>]+W[J] = fj2 J »(x)-Ua(x) (2. 
J* a 

W[J] = lnZ[J] . 



20) 
(2.21) 



Strictly speaking the definition is the convex envelope T[u] — 
minj(J x J2 a Ja{x) ■ u a (x) — W[J]). Here we apply the defi- 
nition to the replicated action, and will content ourselves with 
the differential definition 



ST[u] 

Su a (x) 

SW[J] 

SJ a (x) 



Ja(x) 

u a (x) 



(2.22) 
(2.23) 



which relates a pair of values (J,u), later also denoted by 
(J[u),u). Since T[u] defines the renormalized vertices, its 
zero momentum limit defines the renormalized disorder. Thus 
in order to compute the renormalized disorder, we only need 
to compute T[u] (per unit volume) for a uniform configu- 
ration of the replica field u a (x) = u a = VlVv a (a so- 
called fixed background confi gura tion). Because of the sta- 
tistical tilt symmetry I5ll l52l l53tl . i.e. invariance of disor- 
der term in the replicated action J2.16> under the translation 
v a (x) — > v a (x) + w(x), and of the O(N) invariance one can 
argue, and this is what we find below, that for the model ( 12. 41 
the scaled effective action per unit volume (which for a uni- 
form mode is simply a function of u a ) should have the follow- 
ing form 



(2.24) 



where L d is the volume of the system, and here and below we 
use the notation: 



vv := v a ■ v b 



(2.25) 



for the n by n replica matrix. This defines the renormalized 
disorder. Furthermore, whenever U(vv) can be expanded, up 
to a constant, in the form: 



^(H = W~2 E Z{vlb) - i E ^ Via) - 



2T 2 >> abJ 3 | T 3 

ab abc 



where here and in the following we denote 



Vab '■= V a - Vb 



(2.26) 



(2.27) 



then ( I2.26l > defines the renormalized cumulant functions B(z) 
etc.. As we will see below this is correct up to some very 
subtle behavior at coinciding replica vectors (i.e. v a b = for 
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some pair a, b). Also note that the constant part U(v.v = 0) 
is the free energy. 

The main result of the following sections will be the exact 
calculation of the uniform part of the effective action, i.e. of 
the function U(vv). This will be performed within a large N 
expansion: 



U(vv) = U (vv) + —Ui(vv) + 



(2.28) 



and here we will obtain the dominant order Uq(vv); the cor- 
rections Ui(vv) are calculated in i47ll . It will be a function of 
a scale parameter. We choose to add a mass-term m which 
provides such a scale. It is a convenient choice since for 
rn = oo one has U = U: Fluctuations are totally suppressed 
and the effective action equals the action. One can then pro- 
gressively lower the mass down to zero, starting from this ini- 
tial condition, since ultimately one is interested in the mass- 
less limit. Another choice is to change the UV cutoff, as will 
be discussed again below. 

It is now useful to give a more direct physical interpreta- 
tion of this quantity, in addition to the above field theoretic 
interpretation. 

2. Effective action as the distribution of the order 
parameter 

The effective action for a uniform background is also known 
to be related to the distribution of the order parameter. Let 
us recall the relation for a simple pure ferromagnet. The 
unnormalized probability distribution of the order parameter 
$ = -p- J <f)(x) where <fi(x) is the local magnetization is by 
definition 



Z($) = / V[<f>] 6 $ 



L" 



(x)\e 



(2.29) 



where S[tf>] is the action which describes the ferromagnet (e.g. 
a <j) A theory or a Landau Ginsburg model). The functional 
W( J) evaluated for a uniform J reads: 



W{J) = J d$Z($)e id 



d$ e ^ 



(2.30) 



In the large-volume limit, the saddle point can be taken and 
since the Legendre transform is involutive, this yields the re- 
lation between the effective action at q = per unit volume 
and the probability distribution of the order parameter as: 



-f($)= lim ^lnZ($) 



(2.31) 



In the thermodynamic limit the effective action per unit vol- 
ume can very well be a non-analytic function. This is the case 
e.g. in the ferromagnetic phase where its left and right second 
derivatives at $ = M do not coincide (M is the spontaneous 
magnetization per unit volume). While the right derivative at 
<f> = M is related to the inverse susceptibility, the left one is 



zero, mathematically due to the prescription to take the con- 
vex envelope, and physically because one can always lower 
the magnetization at no cost in free energy per unit volume 
by introducing a domain wall. The above property J2.3H can 
be extended to a given q mode. Finally, note that in d = 
the above does not hold since there is no large factor L d , 
and the probability distribution is directly given by the action 
5(* = <j>). 

What is then the physical meaning of the quantity that we 
will be computing in the next sections? Let us in analogy to 
the magnetization for a ferromagnet define the center of mass 
of an interface: 



1 



N l / 2 L d 



u{x) . 



(2.32) 



Since we have added a mass in the elastic energy ( 12. 4> . which 
acts as an extra quadratic well bounding the fluctuations of 
the interface, the disorder-induced fluctuations of the center of 
mass are always finite. One expects that they diverge typically 
as w ~ 77i~£ as m — > 0, thus their behavior as a function of 
m is of high interest and yields e.g. the information about the 
roughness exponent. 

One can then define the probability distribution Pv{w) of 
the center of mass of the interface in a given realization of 
the random potential V (and in presence of the quadratic well 
induced by the mass). One can see that by definition the gen- 
erating function for a uniform j is the Laplace transform of 
the probability distribution of w, namely 

Z(j) = J dvn . ..dw n P v ( Wl ) . . . P v {w n )e- NLd ^«^ w « , 

(2.33) 

then by the same saddle point argument as for the ferromagnet 
one expects, at least naively, that 

HW] = - Hm^ -Lj In P v (w 1 )...P v (w n ) . (2.34) 
Symbolically one can write: 

Pv(w 1 )...Py(w n ) a e - LdNt ^ , (2.35) 

provided this is taken with a grain of salt. Thus one can also 
think of the renormalized disorder NU(v ■ v) as parameteriz- 
ing the set of correlations of an effective equivalent toy model 
( d = 0) which has the same set of correlations as the center of 
mass variable in the original model. 

The p-th connected moment of the center of mass is identi- 
cal, up to a volume factor, to the zero momentum limit of the 
connected m point correlator of the u field, e.g. 

(w ai ...W ap ) c = ^ (vaAll) ■ ■ ■ v a p {q P )) c \q z =Q (2.36) 

and, once the effective action is known, both can thus be ob- 
tained in principle as the sum of all tree graphs made from 
T[v) vertices. For instance the 2-point function should be ob- 
tainable from: 



{WaWb) 



^- d G ab {q = Q) = ±- d [f"[v = Q]]-J (2.37) 
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and the connected 4-point function from: 

(w a w b w c w d ) con = j^G abcd (qi = 0) 



= ^E[r"b = o]]ae 1 [f"b = Q]] 6 7 

efgh 

x[T"[v = 0}]-i[r"[v = 0}]^ 
xf""[v = 0] efgh (2.38) 

this however assumes analyticity, which as we will see below, 
does not always hold. Another integral relation holds 



(w a w b ) = J dwi ■ ■ ■ dw n wiw 2 Pv{wi) . . . Pv{w n ) 

dwi ■ ■ ■ dw n wiw 2 e~ LdNt [{Wa}] . (2.39) 



III. CALCULATION OF THE EFFECTIVE 
ACTION 

Let us now consider explicitly the large N limit. One can 
rewrite for any N the starting generating function ( 12. 1512.161 

as: 



Z[J] = J V[u\D[ X \D[\]e- NS ^ x ^ 



(3.1) 



1 



<%,x,A,j] = -J C(q) 1 v a (~q) ■ v a (q) -J j a (x) ■ V a (x) 

+ J u (x(x)) - ^iKb(x) [xab{x) - v a {x) ■ v b (x)] , (3.2) 

where the replica matrix field x( x ) = Xab{x) has been in- 
troduced through a Lagrange multiplier matrix \ ab {x). Here 
and below summations over repeated replica indices are im- 
plicit. One can then explicitly perform the functional integra- 
tion over the field u(x) and obtain: 



Z[J]= J D X D\e 



-NS[ x ,X,j] 



(3.3) 



S[x,X,j] = -Trln(C- 1 +i\) 



+ I U(x(x)) -V\x) X ab {x) 



j a {x){C- l +i\)-l bxl j b {x'), (3.4) 



where the inversion and trace are performed in both replica 
space and spatial coordinate space. 

It has now the standard form for a saddle point evaluation of 
the functional W(J) = lnZ(J) except that the saddle point 
is not, in general, uniform in space. It is useful to define the 
scaled functional W(j) through 



(3.5) 



W[J] = NW[j = J/VN] 



which has a well defined large-TV limit and can be expanded 
in 1/N as: 



W\j} = W°\j] + ±W 1 \j] 



(3.6) 



Deferring the calculation of the corrections to a future pub- 
lication 14711 . we obtain here the dominant order in 1/iV as: 



W°[j} 



(3.7) 



where Xj an d depend on j(x) and are the solutions of the 
saddle point equations obtained respectively by setting to zero 
the functional derivatives (at fixed j(x)): 



SX ab (x) 
5S[x,\,j] 







(3.8) 
(3.9) 



X=X:i 



,A=A, 



The result is 



xf(x) = {Gj) ax , bx + {Gj : j) ax ■ {Gj : j) bx (3.10) 



i\f(x) = 2d ab U( X j(x)) 



C- 1 +iX, 



(3.11) 
(3.12) 



where Gj is a matrix with both replica indices and spatial co- 
ordinates and inversion is carried out for both. Here and be- 
low, replica indices are raised whenever explicit dependency 
is given, e.g. Xab = x'f ■ The notation for the iV-component 
vector (G : j)\ x = J2 c J y G bXiCy f c (y) is a shorthand for a 
matrix product, and everywhere we denote by 



d ab U((f>) := dfaUtf) 



(3.13) 



the simple derivative of the function U((f>) with respect to its 
matrix argument (f> ab . Of course, if, for a given j(x) there are 
several solutions to these equations, then one must sum over 
all saddle points, to the same order 



W[j] « In 



sp 



-Ars[xs P (i),A 8p (i),i] 



(3.14) 



This case will be discussed below, for now we ignore this pos- 
sible complication, as well as issues of stability of the saddle 
point. 

Now we want to take the Legendre transform and trade the 
variable j for the variable v to obtain the effective action T [u] . 
One also defines the scaled functional, and its 1/7V expansion 
through: 

~ (3.15) 
(3.16) 



I>] = NT[v = u/VN] 
r[v]=T [v] + ±T 1 [v] 



Then (f [w],VK[j]) and (f o[w],Wo[j]) are also two pairs of Leg- 
endre transforms. Thus the dominant order of the effective 
action functional in the large- N limit is given by 



r [v] = \ v a (x)-jS(x)-w [j v } 



(3.17) 
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with Wo [j] given by (13. 71 . (13. 121 and where (x) is the v(x) 
dependent source solution of: 



5W [ Jv ] 
5j%(x) 



= v a (x) 



(3.18) 



One can now derive a self-consistent functional saddle point 
equation for Tq [v] . First we establish the relation between v 
and j v , namely 

v a (x) = (G v : j v ) ax f v {x) = (G' 1 : v) ax , (3.19) 

where from now on we define 

G v := G jv . (3.20) 

Eq. d3 . 1 91 is obtained noting that 



v a (x) 



SWo\jv\ 



d 



dxM dS | dXjjy) dS 



dj a (x) dxj(y) dj a (x) d\j(y) 



OS 



dja(x) 



dS 



dja(x) 



{Gjv '■ 3v)c 



(3.21) 



where we have used the saddle-point equations ( 13. 8> . ( 13. 9> . 

We can now use ( 13 . 1 9I > in the saddle point equations ( 13. 8L 
(13.91 and defining 



A) 



A; 



Xv '■ Xjv 1 ' y V ■— "»J„ ? 

this yields a self-consistent equation for Xv(x) 

Xv b (x) = V a {x) ■ Vb(x) + (G v )ax.bx 



(3.22) 



(G^ 1 )ax,by = {C-%, y S ab + 2d ab U{ X v{x))S d (x - V ) , 

(3.23) 

which is also a self-consistent equation for G v . Since the Leg- 
endre transform is involutive, one can also write: 



sr Q [v] 

Sv a (x) 



j a {x) = [G v 1 : v) ax , 



(3.24) 



which determines the derivative of Tq [v] once J3.23I is solved. 

One can however do better. Using d3 . 1 91 in ( I3.17> one ob- 
tains the effective action for a spatially varying field v(x): 



T [v] = v : (G v 1 ) : v + S[x v , K, jv] , 
which gives 



(3.25) 



r M 



X1J 



G ax,by V a{x)v b {y) 



+ ±TvMC- 1 + 2dU( Xv )) + j U( Xv (x)) 

+ / v a (x)d ab U(xv(x))v b (x)-xf(x)d ab U(xv{x)) ■ 

J X 

(3.26) 



It is interesting to rewrite it with the help of d3.23i as a func- 
tional of G v and v only: 

T Q [G v ,v] :=T [v]=-^TrlnG v 



+ 9/ Cj by [v a {x)v b (y) + (G v ) ax ,by 
L Jxy 

+ / U(vv(x) + (G v ) x . x ) . 



(3.27) 

We have dropped a constant ~ n. (I3.27> has the property 

d G r [G,v}\ G ^ =0, (3.28) 

where the derivative acts only on G, leaving fixed all v, since 
this coincides with the saddle point equation ( I3.23> . This 
makes apparent that it can also be obtained from a vari- 
ational method where the average of the field is fixed, as 
we detail in Appendix IXI Since the explicit non-trivial v- 
dependence in d3.27i using d3.23i is purely in terms of the 
bilinears v a (x) ■ v b (x) at the same space points, it also shows 
that one can write: 

r [v] = ^v:C- 1 :v + U [vv], (3.29) 

where the interaction (i.e. disorder) part satisfies 

— — -— = , x y (3.30) 

d(Va{x) ■ Vb{y)) 

and is the solution of a self-consistent functional equation: 

SUq[v ■ v] 



6(v a {x) ■ V b (x)) 



d ab U(vv(x) + (G v ) x . x ) 



(3.31) 



(r<-^\ - (n-i\ x SUq[v ■ v] d 

' )axM ~ (C )x - ydab + 5(v a (x) ■ v b (x)) " [X _ V} ■ 

A generalization of this equation is presented in AppendixlBl 

IV. SELF-CONSISTENT EQUATION FOR 
THE RENORMALIZED DISORDER 

A. Uniform configuration and saddle point 
equation 

Let us now consider the simpler problem of computing the ef- 
fective action for a uniform field configuration, which can be 
solved self-consistently. To be more specific we will focus on 
the form ( 12. 4t for the elastic energy. Also, to simplify nota- 
tions and since we will restrict ourselves to dominant order in 
X/N, we drop the index 0, so we set: 



r ^r , U ^U 



(4.1) 



and so on. 

For a uniform field v a (x) — v a the effective action ( 13.291 1 
per unit volume takes the form: 

f(v) := Jjf (v) = ^m 2 v 2 a + U(vv) . (4.2) 
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Note that these are now simply functions (not functionals) of 
a N x ri-component vector, and U(vv) is a function of the n 
by n matrix v a v b . 

Eq. J3.27i yields also a formula for U~o(vv) (up to a con- 
stant): 

U(vv) = U( X v) + \J q * I { H(l 2 + m2 ) S + 2TdU( X v)) 

+ (q 2 + m 2 ) [(q 2 + m 2 )8 + 2TdU{ Xv )] ~* } .(4.3) 

The trace acts in replica space, and the log is a function of a 
matrix, to be defined as usual. Since (14.3b contains the deriva- 
tive d ab U°(vv) we must first determine the latter. One finds 
that analogous to (I3.3H 



d ab U(vv) = d ab U{xv) 



(4.4) 



X f = v a v b + T / [{q 2 + m 2 )5 + 2TdU( X v)]- b L (4.5) 

Jq 

Since one can replace the matrix dU(xv) by dU(vv) in the 
denominator of ( 14.51 1. this is also a self-consistent equation, 
which involves only d ab U(vv). Here inversion is simple n 
by n matrix inversion and S is the Kronecker n by n identity 
matrix 5 ab . One must be careful that 



d ab U{vv) 



8U(vv) 

d(v a ■ Vb) 



(4.6) 



is a first derivative of U(vv) with respect to the matrix ele- 
ment v a ■ v b . One can also check that taking the derivative 
of ( 14. 3> with respect to v a ■ v b correctly reproduces ( 14. 5> . A 
direct derivation uses dxd(vv) from Eq. ( 14. 5> . A more clever 
way is to remember that because of (I3.28> . one is allowed to 
differentiate only with respect to the explicit vv in Xv in the 
first term, and that the remaining terms can be written as a 
function of G v only, and using again (I3.28> . 

This self-consistent equation for dU(vv), i.e. for the uni- 
form part of the effective action is one of our main results and 
the remainder of this paper is devoted to analyze it. It contains 
a huge amount of information, since it encodes the full distri- 
bution (i.e. all cumulants) of the renormalized disorder, and is 
thus quite non-trivial to analyze. It includes both the Gaus- 
sian variational Method (GVM) of Mezard-Parisi|24] and the 
functional renormalization group (FRG). For simplicity, we 
now consider the bare disorder to be gaussian and set all bare 
cumulants except the second cumulant B(z) to zero. 

The GVM is recovered upon setting v = which is one 
limit in which the equation "simplifies". One sees that ( 14. 5> 
then reproduces the Mezard Parisi equations, the self energy 
a ab and two point function G ab (k) in Ref. 12411 being: 



a ab = 2Td ab U(xv=o) 
G ab {k) — G^ = o(k) 

(Xv=o)ab = / G ab {k) . 

Jk 



(4.7) 
(4.8) 

(4.9) 



In the glass phases, these exhibit spontaneous replica sym- 
metry breaking (RSB) with multiple solutions corresponding 



to saddle points obtained via replica permutations, and the 
above equations are solved by a hierarchical Parisi ansatz for 
x(v=0) ab — x( w =0)( u ) where < u < 1 is the overlap be- 
tween replicas a and b. We will give more details about this 
correspondence in the following. 

For now we will study the opposite limit of "strong" explicit 
symmetry breaking field (all v ab = v a — v b ^ 0). Then we 
expect that the renormalized disorder U{vv) is given by a sin- 
gle saddle point and can be expanded in replica sums in terms 
of unambiguous renormalized cumulants, i.e. up to a constant 



1 



3!T 3 



^S(v 2 abl vl cl v 2 ca 



(4.10) 

This is the limit solved here, which we will show below is the 
natural limit in the FRG, and amounts, as we will discuss, to 
forcing the manifold in distant states within the RSB picture. 
The rich crossover to RSB contained in fl4.5ft . when some of 
the v a b are set to zero will be discussed below. 

B. Cumulant expansion 



We now transform equation J4.5I > for the formal function 
U(vv) in a set of equations for the second, third, fourth, a.s.o. 
cumulants. This is performed through an expansion in sums 
over an increasing number of free replica indices, and is not 
an approximation. The such obtained equations are as exact 
as ( 14. 51 . i.e. exact to dominant order at large N, albeit more 
explicit. In fact, they allow a recursive exact calculation of 
all cumulants. Their increasing complexity will illustrate the 
wealth of information summarized in J4.5I >. 
Let us first rewrite J4.5I > using an infinite series: 

dabU(vv) = d ab U( X v) (4.11) 

oo 

Xv = vv + Th5 + T In+i (-2TdU(vv)) n 

71=1 

1 



In ■ = 



k (k 2 + m 2 ) n ' 



(4.12) 



where the n-th power here denotes the matrix product. 

Since we consider a gaussian bare model ( 12.181 1 where only 
the second cumulant is non-zero one has: 



-2Td ab U{ Xv ) = t 



-B'(x a v b ) + SabY, B '(Xv C 



(4.13) 

using that d a bXab = 2(S a b — 1) 15411 . The same quantity for 
the renormalized disorder reads: 

- 2Td ab U(vv) = | {-B' ab + 5 ab ]T B' ac ) 



2 



(4.14) 



where we denote B' ab = B'{v 2 ab ), S abc = S(v 2 ab , v 2 bc , v 2 ac ) and 
S[ abc denotes a derivative with respect to the first argument 
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of the function S (S has the symmetries implied by replica 
permutation symmetry). All matrices we will encounter can 
be parameterized as: 



X a b 


= Xab 4 


Sab%a 




(4.15) 


%ab 


_ r (0) 
— x ab 




r (2) + 

X ab ' ' ' ' 


(4.16) 


■'■a 


- T (0) . 




X& + 

■''a ' • • • I 


(4.17) 



where x ab do not contain any explicit Kronecker 5 ab , the 
upper index denotes the number of free replica sums, e.g. 



X/ f x a 



/ X ab;f, 



,.(2) 



J2f g x ab;fg- Since under matrix 
product (X p ) ab or Hadamar product (X ab ) p the number of 
sums can only increase, one gets only a finite number of terms 
in projecting out on terms with a given number of free replica 
sums. 

If we parameterize in the same way: 



Xv = Xab + SabXa 
Xv Xab $abXa 



(4.18) 
(4.19) 



then one easily sees that: 

B'(xt ) = 5ab[B'{ X aa + Xa) ~ B> ( X aa)] + B' ( X ab) 

B'(xt) ~ Sab £ B'{XT) = B'{ X ab) - Sab £ B> ( X ae) . 
c c 

(4.20) 



We can now expand in number of sums: 

(4.21) 

and the equivalence of (I4.13> and (I4.14> using 14.201 . implies: 



B\Xab)=B\Xa a b ) )+B"(Xa"b ] )Xa 1 b ) 
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B' ab = B'(x<$) 

ibg 



iVs" - B"Cv (0) lv (1) 

-r / , °\,abg ~ D \X a b )X a b 



(4.22) 
(4.23) 



and so on for higher cumulants. Thus to obtain the second 
renormalized cumulant we only need to compute the part y^v 
of Xy b which contains zero sum and no explicit S ab . One has 
in general: 



X a P b ] = + XS } + xi P) + X { b P) - 2x1^ ■ (4.24) 



Thus for the second cumulant we need only x a b an d Xa- Since 
one has, to be explicit: 



J 



[{-2TdUf] a b = 55 



Sab ^ B'aeB'af ~ B' ab ^^{B' af + B' bf ) + ^ B' ac B' ct 



ef 



4 

"T3 



2S ab ^ B' ae S' agh — B' ab ^^{S' bgh +S' agh ) ~ ^^{B' af ,S' abh + B' be S' abh ) + *^^{B' ac S' cbh +B bc S' cah ) 

egh gh eh he 

(4.25) 



where all terms not written have at least three free replica 
sums (this is the case for 0(S 2 ) as well as terms involving the 
fourth cumulant and higher). Similarly [(— 2Td x U) p ] ab has at 
least p — 1 free replica sums (from the 0(B P ) term). This is 
much more what we need, which comes only from J4. 14i and, 
using ( 14.1 1> : 

X ( ab = v a v b - 2I 2 B' ab (4.26) 
Xa 0) = Th . (4.27) 

This yields: 

X^ - (v a - v b ) 2 + 2Th + AI 2 (B' ab - B' aa ) . (4.28) 

Thus we find that the renormalized second cumulant satisfies 
a closed equation at any T: 

B'(v 2 ab ) = B'(v 2 ab + 2Th + 4h(B'(v 2 ab ) - B'(0))) (4.29) 



with no other contributions from higher cumulants at any T. 
Appendix iDl contains a non-local extension of this formula. 
Eq. ( 14.291 can be integrated with the result 

B{v 2 ) = B(v 2 + 2Th + AI 2 [B'{v 2 ) - B'(0)]) 

-2I 2 {B'(y 2 + 2Th + 4I 2 [B'(v 2 ) - B'(O)f) }' . 

(4.30) 

A direct derivation from (13.31 is also possible. 

C. Higher cumulants 

Higher cumulants of the renormalized disorder can be ob- 
tained by the same method using ( 14.231 and its extensions. 
They can also be obtained by the graphical method. For sim- 
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plicity here we give only the expression of the third cumu- 
lant. The complete expression for the fourth cumulant, to- 
gether with all calculational details and an introduction to the 
graphical method, can be found in Appendixlcl 
The third cumulant is found to be: 



S(x,y,z) = 



6TI 2 



: s y m x, 



B'{x)B'{y) 



1 + 4/ 2j B"(0) 
+24/ 3 synv^[(B'(x) - B'(0))B'(x)B'(y) 
-8I 3 B'(y)B'{z)B'(x) , (4.31) 



where sym^, z is 1/6 times the sum of all permutations of 
x,y, z. Note that this relation is exact for all values of the 
mass m, and not just a fixed point form. The only input in 
the derivation is the absence of a third cumulant for the bare 
model (m = oo). It would be interesting to include an addi- 
tional bare third cumulant. The fourth cumulant is derived in 
appendix |C] where also details for the graphical method are 
given. 



V. GRAPHICAL INTERPRETATION 

In this section, we sketch how the central results at large N 
can be obtained graphically, first the saddle-point equation 
d4.29i . which gives the effective disorder B as a function of 
the bare disorder B, and second the /3-function ( 16. 9i . 

The graphical rules for the perturbation theory of the repli- 
cated model have been described in detail in 15511 for N = 1 
and we refer the reader to this work for elementary details. 
Here there are in addition N components of the field u l a , the 
propagator being diagonal in all indices. For the present pur- 
pose we are mostly interested in the counting in N, and since 
it is difficult to represent graphically both vector- and replica- 
indices, we work with unsplitted vertices (see 155]) and spec- 
ify the replica content only when needed. Disorder vertices 
may contain arbitrary number of derivatives and some exam- 
ples are represented on Fig. |2] As usual there is a factor of 
1 /N per derivative (i.e. per dashed line) at each vertex (see 
e.g. Fig. |2 using v 2 = u 2 /N), N per vertex, and N per loop. 

We consider the effective action, i.e. the sum of all 1- 
particle irreducible diagrams (1PI), and later focus on its 2- 
replica part. We start our analysis at T = with the three 
possible 1-loop diagrams, as presented on figure [5] They are 




FIG. 1: Graphical representation of the third cumulant. The notation 
is explained in 1 47] . The first diagram yields the terms proportinal to 
I2, the second diagram the terms proportional to J3 in Eq. J4.3 II . 



; r 



-B (V )(« ) 



±B"'(v 2 )(v 2 ) 3 





I 

TB'{v 2 ) 



FIG. 2: Examples for vertices and the 1-loop tadpole diagram which 
is dominant at large N. 






FIG. 3: The four 1-loop diagrams correcting the disorder. A fat dot 
represents a vertex B, a solid line the field u, and its correlator. A 
dashed line attaches two fields u to a vertex B. We do not draw 
replica-indices. 



obtained from contracting 



N 



t2>( ! 



N 



(5.1) 



ab 



cd 



In order to simplfy the calculation we omit the terms taken at 
coinciding replicas (e.g. B'(0)), they can be added at the end. 
Contracting j5.U twice between points x and y gives 



^c 2 V 

9 xy 



B'( 



i( {<-0 2 \ B n< «-<) 2 \ K-O 2 



N 



ab 

±B 

N 

A B 

N D \ N I N 



N 



N 



)1 B '< ( - 



b\2 



N 



N 



A' 



■B 1, 



N 



N 



higher replica terms 



(5.2) 



This is graphically depicted on figure [3] The important ob- 
servation is that only the first diagram, with a closed u-loop 
is contributing in the limit of large N. This analysis can be 
repeated to higher loop order. Again, only diagrams as the 
first one on figure |3]contribute. Especially, there are no loops 




o 



FIG. 4: Loops which give additional factors of 1/N, as explained in 
the main text. 
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FIG. 5: Tree-configurations which contribute to B(v 2 ). 



U | U i I) 0=b(x) 

2 d 6 o "' • 

o o 

FIG. 6: Self-consistent equation at leading order for B'(v 2 b ) — 
B'(xab)- The wiggly line denotes a derivative, and is combinato- 
rially equivalent to choosing one B. At finite T one can attach an 
additional arbitrary number of tadpoles to any B. 



with three propagators or more, as loops 4 or 6 on figure 0] 
Also, there are no "meta-loops", i.e. loops formed by loops, 
as loop 5 on figure^] Finally, only diagrams as those on figure 
|5] survive, which as building block have only the elementary 

1- loop diagram with a closed loop contributing a factor of N, 
as the first diagram on figure [3] These are tree-like diagrams, 
where the nodes are made out of B and the links out of the 
before-mentioned 1-loop diagram (two parallel replica lines 
in the splitted diagrammatics [55] which produce the desired 

2- replica term). At junction points the replica lines branch 
also in parallel. These are of course not tree diagrams, i.e. 
they are 1PI and contribute to the effective action. Note on 
Fig. |5] that since there is a 1/(2T 2 ) factor per vertex, but that 
each vertex (except one) comes with two propagators (factor 
T 2 ) the counting in temperature is right to produce a 2-replica 
term with the expected 1/(2T 2 ) global factor (the 3-replica 
terms proportional to T have been discarded etc.). 

We are now in a position to derive the self-consistent equa- 
tion ( 14.291 . The key-observation is that deriving B(v 2 ) once 
with respect to its argument, amounts in the graphical in- 
tepreation of figure|5]to choose one of the bare vertices B(v 2 ), 
and deriving it. B'(v 2 ) thus is B'(v 2 ), with as many branches 
attached as one wants. Every branch consists out of a 1-loop 
integral <^> = I 2 times another tree; the latter is again 
B(v 2 ), given that one of the bare vertices is chosen, i.e. again 
B'{v 2 ). Since attaching loops to B' [v 2 ) amounts to deriving 
B'(v 2 ) once for every loop, we arrive at 



oo B^ n+1 \v 2 ) AI 2 {B> \v 2 ) - & '(0)) 

£=0 



B 



'(v 2 +AI 2 [B l {v 2 )-B'(Q)] 



(5.3) 



Note that we have added the term with coinciding replica- 
indices, dropped previously. The combinatorial factor comes 
from the expansion of the exponential function in e~ s . That 
it indeed resums to B' with a shifted argument is natural: For 
a function f(x) taking the expectation value (f(x)) in a the- 
ory with only a first moment (x) is equivalent to calculating 



/((&)). Taylor-expanding the latter leads to the above combi- 
natorics. 

By the same arguments the full effective action can be writ- 
ten as the sum over tree-like (but not tree) diagrams repre- 
sented in Fig.[5]where, in addition, each vertex can be dressed 
by an arbitrary number of tadpoles (see Fig(3J- Each tadpole 
brings an additional factor of T, thus tadpoles contribute to 
the two replica term only at T > 0. At finite temperature, any 
of the v 2 's could be contracted, leading to the relacement 



,,£n2 



„b\2 



2Th 



(5.4) 



(This offers another possibility to verify the combinatorics in 
15.31 ) Thus the final result is 



2Th+M 2 [B'{v 2 )- B'(0)}) (5.5) 



B'(v 2 ) = B' 



We now illustrate how to recover the /3-function. Applying 
—md/dm to B implies to derive each integral w.r.t. m ap- 
pearing in each loop of Fig|5] Diagrammatically this amounts 
to choosing in the tree of Fig |5] one of the bonds (loop I 2 ) 
which connects two B's. Summing over all trees, it gives a 
term 



9 J 
-m—h 
om 



B'{v 2 ) 2 - 2B'(v 2 )B'(0) 



(5.6) 



since the two trees attached to the loop I 2 are nothing but 
B(v 2 ), derived once, and again itself with things attached, 
i.e. B'(v 2 ) as given in ( 15 .51 . This reproduces the T — term 
in 

The second contribution comes from deriving TI\. The 
graphical derivation is complicated, and we refer the inter- 
ested reader to l47ll where a more complete, but much more 
involved, diagrammatic method is presented. 

VI. FUNCTIONAL RENORMALIZATION 
GROUP EQUATIONS 

A. From self-consistent to FRG equation 

We will now study the self-consistent equation, exact for N = 
oo, for the second cumulant correlator of the random potential 
that we have derived in the previous Section: 



B\x) = B'[x + 2Th + AI 2 {B'{x) - B'(0)) ) , (6.1) 
which involves only the two one loop integrals: 



h 



h = 



A d d k 1 
o {2n) d k 2 + m 2 

A d d k 1 
o (2ir) d (k 2 + m 2 ) 2 1 



(6.2) 
(6.3) 



where we have indicated symbolically that a short scale UV 
cutoff is needed for I 2 to be finite if d > 4 and for 1\ for 
d> 2. 

There is a simple way to obtain directly the solutions of 
( 16. 1> which we will detail below. It is also interesting to turn 
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this equation into a FRG equation for the function B(x) as a 
function of the scale parameter m. Indeed this yields the (3- 
function of the field theory in the limit of infinite N, which is 
our main goal. Let us show first how one does it. 

Let us first take a derivative of i6.1\ with respect to x. One 
obtains: 



B"[x) 



= B" x 



2TI 1 +4I 2 (B'(x)-B'(0))j . 

(6.4) 



l + 4I 2 B"(x) 

Taking the derivative md m of J6. It and using ( 16. 41 gives: 

B"{x) 



md m B'(x) 



2md m T I\ 



l + 4I 2 B"{x) 
+4(md m I 2 )(B'(x)-B'(0)) 
+4I 2 md m B'{x) - 4I 2 md m B'(0) 

Regrouping the terms one obtains: 

md m B'(x) = B"(x) \2md m Th - 4I 2 md m B'(Q) 
+ 4{md m I 2 ){B'(x)-B'(0)) 

From ( 16. 5> one also has 

B"(0) 



(6.5) 



md m B'(0) 



1 + 4I 2 B"(0) 
Inserting ( 16. 7> into (16.61 ) finally yields 

md m B'{x) = B"(x) \2(md m Th) 



2md m {Th 



.(6.6) 



(6.7) 



1 



l + 4/ 2 B"(0) 
+ 4(mdmh)(B'(x)- B'(0)) 



.(6.8) 



This equation is valid for any space dimension d. It can be 
integrated once w.r.t. x to obtain the final result 



mo rn B(x) = = B (x) 

V ' 1 + 4I 2 B"(0) 

+2(md m I 2 ) \B'{xf - 2B'(0)B'(x) 



(6.9) 



where we have dropped an m-dependent integration constant. 

A general method to study (and solve) the FRG equation 
J6.8I is then to start from to = oo where the initial condition 
is B(x) — B(x) in the presence of a UV momentum cutoff A, 
or a lattice with lattice constant a = 1/A. Then one studies 
how B(x) evolves as m is slowly decreased. 

There are thus two possible paths to solve the problem, 
namely the direct inversion of the self-consistent equation and 
the solution of d6 . 8t > with the above initial condition. Both 
are studied below. These two methods are clearly equivalent 
when the solution B(x) is analytic at x = 0. Indeed, in the 
above derivation, we have assumed that B"(Q) exists. This 
will not always hold, as we now discuss. What the proper 
ensuing modifications are is a subtle point which will be ex- 
amined later. 



B. General features: Analytic vs non-analytic 
solution 

Before solving this equation let us first find the conditions un- 
der which there exists an analytic solution. This will give us 
insight in the phases of the model. One notes from ( 16. 41 that: 



1 



B"(0) B"{2Th) 



4h 



(6.10) 



For m = oo the starting value is B"{{)) = B"(0) > 0, in 
any dimension d. (The force correlator decays for small dis- 
tances.) As to is decreased several things can happen. 

Let us start with T = 0. Then for d < 4, since I 2 diverges 
for small to, one sees from J6. 10b that B" (0) becomes infinite 
as to — > to+, where the Larkin mass m c is the solution of: 



„d-l 



iS d / dq 



1 



(<? 2 +to2) 2 B"(0) 



(6.11) 



with S D = l/(2 d " 1 7r d / 2 r[d/2]) and has the standard depen- 
dence m c - B"^) 1 ^ of the inverse Larkin length on the 
bare disorder (a Larkin length L c = l/m c can be defined). 
Since B"(0) is like R""(0) positive, this divergence is the 
usual one of the FRG, as also found in 1- and 2-loop stud- 
ies 1351 l36t l43l 144. 45, 56], where it signals that the function 
R(u) becomes non-analytic and that a cusp singularity forms 
at u = in the second derivative —R"(u), i.e. in the correla- 
tor of the pinning force. This is usually interpreted as a glass 
phase with many metastable states beyond the Larkin length. 
Thus for d < 4 the function always becomes non-analytic at 
large scale (small mass), and there is a single glass phase. For 
d > 4, since I 2 is convergent, the cusp occurs only if the bare 
disorder is sufficiently large. 

At non-zero temperature T > J6. 10b shows that for 2 < 
d < 4 thermal fluctuations do not change the scenario. Since 
I\ remains finite, temperature only slightly renormalizes the 
value of to c downward, as 



4S d / dq 



1 



o y (g 2 +TO 2 ) 2 B"{2TS d K d - 2 /{d- 2)) 



(6.12) 



for A to c . For d < 2 the effect of thermal fluctuations 
is more important. For definiteness let us consider the set of 
models with power law correlations ( 12.10k Then d6 . 1 01 be- 
comes: 



' —(a 2 + 2Th) 1+ ~i -4I 2 . 



B"(0) gi 



(6.13) 



Since both integrals diverge for small mass as I\ ~ 1/ to 2 d , 
I 2 ~ l/TO 4 ~ d , one can distinguish three cases: 

(i) If disorder correlations decay fast enough 7 > j c (d) = 
2/(2 — d) then the I\ term wins and as m — > one has 
B"(0) — » 0, indicating that disorder is subdominant, 
resulting in a high-temperature phase. In that case the 
solution is analytic as m — > 0. There is however a more 
complicated behavior for intermediate values of m (see 
AppendixlEli. 
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(ii) If disorder correlations decay slower, i.e. 7 < j c (d), 
the term proportional to I 2 wins and the solution always 
becomes non-analytic at some Larkin mass. 

(iii) In the marginal case, 7 = j c (d) there is a transition at 
some critical temperature T c between a high tempera- 
ture phase and a glass phase. 

These features are very general and each of these cases will 
be studied in more details below. 

One can immediately see that the existence of an analytic 
solution for B(u) is in one to one correspondence to the exis- 
tence of a locally stable replica symmetric solution of the MP 
equations. Indeed the condition for the stability of the RS sad- 
dle point is precisely that the replicon eigenvalue be positive, 
namely that 10: 



A re p(p) = l-4J 2 (p)B"(2TJ 1 ) 



(6.14) 



J 2 (p)= (k 2 +m 2 )- 1 ((k+p) 2 +TT1 2 )- 1 (6.15) 

Jk 

be positive for all p. The RSB instability occurs when the 
lowest eigenvalue, which corresponds to p = 0, vanishes. The 
condition A rop (p = 0) = is equivalent to the vanishing of 
(16. 101 . i.e. of the divergence of 73" (0) and the emergence of 
non-analytic behavior. Thus the generation of a cusp in the 
FRG coincides at large N exactly with the instability of the 
RS solution. 

It is easy to see that an analytic solution B(x) of J6. It and 
(16.81 cannot describe the glass phase at T — 0. Indeed when 
B(x) is analytic, Eq. J4.3 II and similar results for higher cu- 
mulants indicate that the full effective action is analytic. It is 
then immediate to obtain correlations from its derivatives. For 
instance, from (I2.36> the 2-point function at q = is simply: 



1 



(u a (q) ■ u b (q)) I g=0 



, r "-uvi// \q— u 9 



T B'(0) 

Oab - 2 7- 



On the other hand, setting x = in ( 16. Il l one finds: 

73' (0) = B'(2Th) . 



(6.16) 



(6.17) 



C. FRG equation for rescaled disorder, d < 4 

The equation d6.8t is valid (for TV = 00) in any spatial dimen- 
sion d. Since one has the exact relation: 



1 



-md m Ii = m I2 



(6.18) 



one sees that the FRG equation J6.8t has a well defined limit 
A — > 00 for d < 4. It makes formulae somewhat simpler so 
we will start by considering this case; the case d > 4 will 
be studied later. Note that although the equation has a well- 
defined limit, its solution may require a UV cutoff (e.g. as is 
manifestly the case in integrating ( I6.18> above). 

Thus from now on we study d < 4 and consider the infinite 
UV cutoff limit. Then one has 



h = A d - 



A d = 



T(Z- d -) 



(6.19) 



with e = 4 — d. It is convenient to define the rescaled dimen- 
sionless function: 



b(x) = AA d m iC -- e B(xm- 2C ) 



(6.20) 



where £ is a fixed number, but for now arbitrary. Note that 
whether one works with B or the rescaled b(x) does not make 
any difference for the possibility of a non-analyticity or a di- 
vergence of the second derivative. 

Then b(x) satisfies the FRG equation in the infinite-N limit: 

- md m b(x) = 0[b] 

= (e-4C)b(ar) + 2C,xb'{x) 



+ h>(x) 2 -b>(x)b'(0)+T m - 



6"(0) 



(6.21) 



The rescaled temperature, and the energy exponent 8 are de- 
fined as 



T m = 1 m 

e 

9 = d-2 + 2( 



(6.22) 
(6.23) 



Thus at T = one recovers the dimensional reduction (DR) 



result 



-d-2( 



with C = Cdr = (4 - d)/2 instead 
of a non-trivial value for ( expected in the glass phase. Fur- 
thermore since the effective action is analytic, all higher con- 
nected cumulants will trivially vanish at T = (or be equal 
to the bare ones if the bare model contains such higher cu- 
mulants) from the DR property. Clearly, in the glass phase, 
the DR scaling is expected to be incorrect and a non-analytic 
solution should be found, as well as a way to escape ( I6.17> . 
Below we find how such a mechanism occurs within the FRG. 

It will emerge from our study that for the case where dis- 
order is relevant in the large scale limit (i.e. the long range 
case 7 < 7 C (d) mentioned above) the non-analytic solution of 
the FRG equation will correspond to the full replica symme- 
try breaking solution of MR The situation for the short range 
case is more delicate. Both are discussed below. 



To obtain J6.2 li we have also integrated d6.8i once, so there 
is a priori a m-dependent integration constant. 

We emphasize that this FRG equation ( 16. 2li that we have 
derived is valid, to dominant order in l/N, in any dimension 
d < 4 and at any temperature T. In a previous study 1 3 6| Ba- 
lents and Fisher studied another limit: arbitrary N but only to 
first order in e = 4 — d and T = 0. If we consider the dom- 
inant order in TV of their equation, we find that it is identical 
to the T = part of (I6.2H (up to some changes in notation). 
Equation ( 16. 211 however is valid to all orders in e, an impor- 
tant point which the method used in |36] could not address. 
Comparison of ( 16. 211 to our recent 2-loop, i.e. 0(e 2 ) studies 
requires to expand to next order in l/N and is performed in 



requi 



Furthermore ( 16. 2U includes the effect of temperature to all 
orders in e. Expanding the term proportional to T to lowest 
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order in disorder b, one finds the term T m b'(x). This is the 
large- N limit of the tadpole term obtained in the 1-loop FRG 

atr > o EHE3,ESEai: 

jv 

diR(u) = T^2d?R(u) 

i=l 

-> diBiv 2 ) = TB'(v 2 ) + ^v 2 B"(v 2 ) , (6.24) 

where for infinite N the last term drops out. (It appears how- 
ever to next order in 1 /N |33l-) 

The form and the effect of the temperature term in (16.2 1 1 
to all orders in e is radically different from its 1-loop trun- 
cation. Indeed, in the 1-loop FRG the temperature is known 
to smoothen the cusp and render the function R(u) analytic 
in a boundary layer u ~ T m (e.g. for N = 1 
with R""(0) ~ 1/T m . Here however, as further analysis 
confirms below, for 8 > the divergence of 6"(0) is self- 
reinforcing since it kills the term proportional to T m . We find 
that it usually occurs at a finite (Larkin) scale. In the marginal 
case 9 = 0, we will find non-trivial analytic finite-temperature 
fixed points. 

VII. DETAILED ANALYSIS OF THE FRG 
EQUATIONS 

A. Inversion of self-consistent equation 

Let us now show how one can invert the self-consistent equa- 
tion J6.ll>. We first rewrite it in terms of the rescaled correlator 



b'(x) = 4A d m 2<: - e 

xB'(m- 2C (i + -(b'{x) - b'(0)) + 2TW C )) , 

(7.1) 

where in the term proportional to temperature, for d > 2 we 
mean lirriA^oo TI\ choosing a bare temperature T ~ A( 2 ~ d ) 
(this choice is known to be necessary to give a universal and 
finite (3 function, see e.g. the discussion in Ref. i47ln . One can 
of course keep an explicit A dependence everywhere, but that 
leads to needless complications without changing the result. 
The above equation (17. 1> is easily inverted into 



where we define 



y 



4A d m 2 (- e 



+ ~(v - Vo) - T„ 



y = y{x) = -b'{x) 

y = -b'(0) = -AA d m 2 ^B'{Q) 
f m = 2TI im 2C 



(7.2) 



(7.3) 
(7.4) 
(7.5) 



with T m — T m / (2—d) for d < 2, and $ is the inverse function 
of -B' (x) i.e. 



This means in turn that the FRG equation (I6.2H is fully in- 
tegrable, a feature not immediately obvious if one does not 
know that it originates from a self-consistent equation (an ob- 
servation not made in Ref. 1 36]). To better understand this in- 
tegrability property let us show that (I6.2H can be transformed 
into a linear equation. Let us first take a derivative of (I6.2H 
and express it in terms of the new function y(x) i7.3\ 



-md m y = (e-20y+2<:xy'-y'(y-y )+T m ^- r , (7.7) 

1 Vq 

e 

where we denote y' = y'(Q). Converting this into an equation 
for the inverse function x(y) one finds: 

T ex' 

md m x = (e - 2()yx / + 2(x - (y - y ) + ^ (7.8) 

ex Q - 1 

with x'q — x'(yo) f72Tl . We have used that md m y(x) — 
— y' (x)md m x(y(x)) and have canceled a factor of on 
both sides. (The validity near x = beyond the Larkin length 
is reexamined below). 

One recovers now that the general solution of this linear 
equation is ( 17. 1> since it is the sum of the general solution of 
the homogeneous part 



— m 2 ^(j)[ym e 2 ^ 



(7.9) 



where (f> is an arbitrary function, and of a particular solution 



1 



x=-(y- y ) - T„ 



(7.10) 



The y dependence obviously satisfies (17. 8t and for the con- 
stant part to work we use: 



md m y = (e - 2C)y + T n 



(7.11) 



I _ % 

e 

-md m f m = -2(f m + T m . (7.12) 



(-£?')(<%)) = y ■ 



(7.6) 



The first line comes from evaluating ( 17. It at x = and as- 
suming analyticity, i.e. that \im x ^ y'(x)(y(x) — yo) = 0, 
and equality which will not work beyond the Larkin length 
(to < m c ), as found below. 

Now that we have clarified the connections between the two 
approaches (self-consistent equation and FRG) we can try to 
find solutions valid in the small mass limit. To analyze the so- 
lutions of the large- N FRG equation (I6.2H . two approaches 
are legitimate, corresponding to different points of view. The 
first, natural in mean field, is exact integration. But then one 
discovers that the solution becomes non-analytic upon reach- 
ing the Larkin mass. It thus raises the non-trivial question on 
how to continue this solution beyond the Larkin length. Be- 
fore doing so, we will first examine a second point of view, 
more familiar from standard RG arguments. 

B. The FRG point of view: Search for fixed points 

The standard RG approach amounts to construct and compute 
the /3-function of the theory, and then search for a fixed point 
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(function) which describes the large scale physics. Usually, 
finding the basin of attraction of the fixed point, or relating 
arbitrary initial conditions to the final approach of the fixed 
point is an unmanageably difficult task. It is fortunately also 
besides the goal of the RG which is to compute universal large 
scale physics independently of the irrelevant details of the bare 
model. Here, however, because of the large- N limit, we can 
integrate the RG flow exactly and in principle "solve" any bare 
model. Let us temporarily ignore this integrability feature and 
focus on finding the zeroes of the /3-function. 

The /3-function was derived previously within an e expan- 
sion and non-analytic fixed points were found to one loop 
||a,|35lp6] and also to two loops I43l i45ll55ll . In the latter case 
additional "anomalous" terms are present in the /3-function for 
the non-analytic theory to be renormalizable and a meaningful 
fixed point to exist. Viewing the right hand side of (16.2 1 i as 
the large-TV limit of the true /3-function, let us follow the same 
strategy and ask whether we can find non-trivial fixed points. 

Let us study T = and use the equivalent linear form of 
the FRG equation. We want to find the solutions y(x) of 



(e - 2Qyx' + 2(x - (y - y ) = 



(7.13) 



t/o is a fixed number (we want to impose yo = 2/(0)), since we 
are looking for a fixed point function. Keeping yo arbitrary, 
one first tries a linear solution x = ay + b which yields a = 
1/e and b = -y /(2(). Writing x{y) = (y/e) - yo/(2Q + 
4>(y) one finds a homogeneous equation for <p and thus 



/ s y vo 



ay 



2C 

E-2< 



(7.14) 



Imposing now yo = 2/(0), i.e. x(yo) = 0, fixes the value 
of a and one finds the family of zero temperature fixed point 
functions, parameterized by f : 



2/o f 

2C 2Ce 



Vo y ^ 



(7.15) 



Since x > 0, yo > one must have 



< C < 



> and thus 



(7.16) 



The case C — f corresponds to a Larkin random force model. 
For the same reason, we must exclude the branch y > yo and 
thus x* (y) is given by the unique solution of J7.151 with x > 
and < y < yo- Finally, for ( = 0we find the fixed point: 



x = x*(y) = -(y - y Q - y ln(y/y )) 
e 



(7.17) 



An important observation is that all of these fixed points ex- 
hibit automatically the expected cusp. Indeed one finds that 
x'(yo) — 0, i.e. x(y) in (I7.15> vanishes and has a minimum at 

2/ = yo- 



This gives 



1 



2(e - 2C)2/o 



(y - yo) 2 + 0((y - y Q ) 3 ) . (7.18) 



b'{x) = b'{0) + Ay/x~ + O(x) 



with A = J 2(e - 2()\V (0)\, implying that the second 
derivative diverges as x — > + 



b"(x) 



A 



O(x ) 



(7.20) 



Recalling that y = —b'(x) we see that all fixed points with 
C > correspond to a power-law long-range correlator b(x), 
while ( = corresponds to a gaussian short range disorder. 
If we follow the standard RG arguments, we can now sort the 
models J2.6I > into these universality classes. Since for the bare 
model 



B'(z) 



(7.21) 



and since the decay of R(u) in i2.6\ at large u can be argued 
to be identical for B and B (for LR fixed points) we find 

A-d 



C = C(7) - 



2(1+7) 



(7.22) 



(7.19) 



or ( = for short range correlations. These values are valid 
to dominant order in l/N. In O the effect of the 0(1 /N) 
terms in the 1-loop FRG equation was studied, i.e. the correc- 
tions of ( were estimated to order 0(e) and at zero tempera- 
ture. For SR disorder it was found that the result of the GVM 
(i.e. Flory) is corrected by terms ajy exponentially small in N, 
i.e. Csr = C(7 = f + !) + awe + 0(e 2 ). For LR disorder 
with 7 > 7*(iV) the result (I7.22> was found to be uncor- 
rected to 0(e). (The crossover SR to LR occurs at 7* such 
that £(7* ) = Csr)- One can in fact argue that J7.22l > is always 
exact in the LR case (see e.g. discussion in Ref. |55]). 

Several important remarks are in order. First we have found 
the fixed points of the inverted linear form ( 17. 8> of the FRG 
equation. A valid question is whether this is equivalent to 
finding the fixed points of the initial form of the /3-function 
A6.2 II . Second we have found fixed points assuming that 
Tnd m yo = 0. Since this is different from what has been found 
previously in ( 17.1 1> at T = 0, one can ask whether these result 
are compatible. 

These two questions have a common answer. Examining 
more closely what has really been done in this Section, we 
note that it is equivalent to declaring both J6.2 li and 
valid for any x > and interpreting everywhere yo — 2/(0 + ) 
in ( 17.81 and, equivalently b'(0) as b'(0 + ) defined by continuity 
as x — > + . This is legitimate since the transformation from 
(16.2 1 i to d7.8t is certainly valid for x > and we note that 
this answer the second question above since Eq. ( 17. It , i.e. the 
derivative of ( 16. 2U . evaluated at x — > + yields: 

-md m y(0+) = (e-2CM0 + )- lim y'(x)(y(x)-y'(0+)) , 

x->0+ 

(7.23) 

which works both in the regime m > m c where the solution 
is analytic 2/(0 + ) = 2/(0) and in the fixed point regime m — > 
when the cusp has developed and the last term in ( I7.23l > has a 
non-zero limit according to ( I7.19l7.20t . 

We expect these fixed points to be the physically correct 
solutions at small to. We now investigate whether we can 
confirm this by providing the solution at infinite N, for ar- 
bitrary mass to, i.e continue our solution ( 17. H beyond the 
Larkin length. 
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C. Full solution beyond the Larkin mass 

We now show that one can connect the two regimes, i.e. the 
regime for m > m c where an analytic solution exists to the 
asymptotic one, for to — ► 0, studied in the previous Section. 
This can be done here because of the full integrability of the 
infinite- TV limit and provides a rare and non-trivial insight into 
what happens around the Larkin scale. 

It is instructive to start our analysis with the specific power 
law models with LR correlations ( I2.10t . together with the case 
of SR correlations j2.9l l. in the form of a Gaussian. The so- 
lution for an arbitrary bare potential B(z) is more subtle, and 
will be given in Section fVlll Dl and appendixlHl 

For the power law correlators the inverse function in i7.1\ 

is: 



-B'{z) = 



(a 2 + z)i 
For gaussian correlations it is: 



^ z = $ 



M5) 



-l/ 7 



(7.24) 



-B'(z)=ge- Z o z = $(y)=\n(g/y). (7.25) 



We can now insert this result into the general solution i7.2h of 
the self-consistent equation. £ is arbitrary, but the convenient 
choice (to later obtain a fixed point) is £ = £(7) such that the 
to dependence of the first term drops. Let us define 



9 = ^A d g . 
We then obtain, for power law models: 



-l/ 7 



1 



iv-Vo) 



V 

3 J e 
-6'(0) = y = g(m^a 2 + f m )^ 



2C 2 

to 



T 

1 rr 



(7.26) 

(7.27) 
(7.28) 



since we want y(0) = yo i.e. x(yo) = 0. This solution is valid 
for m > m c and the value of b'(0) is the DR result d6.17> . For 
short range disorder the solution for to > m c is 



x = ]xi(gm e /y) + e (y - y Q ) - T„ 



-tf(Q) = y = gm~ 



(7.29) 
(7.30) 



having set ( = in that case. We recall that y(x) = —b'(x). 
Note that the bare disorder is recovered for to — * 00. We have 
kept temperature, but here we discuss only the case where 



0(7) 



d-2 



2C( 7 ) > , 



(7.31) 



i.e. 2 < d < 4, or d < 2 with 7 < j c (d) = 2/(2 - d). In that 
case T m decreases as to decreases, and, as mentioned above 
the role of temperature is minor. 

Let us plot the r.h.s of J7.27K ( 17.291 1 on Fig.0 The curve 
x(y) has the indicated shape in all cases. It cuts the axis x = 
at y = y and has a minimum x'(y c ) = at y = y c with 



Vc 



gl/(l+7)( e / 7 )7/(l+7) 



(7.32) 



independent of to, and y c — e for SR disorder. For to > 
to c the minimum occurs at negative x and the slope at y = 




FIG. 7: The function x(y) given by <T27l or fT29l . The physical 
branch is the one with y < yo- 



yo < Vc is non-zero, indicating an analytic solution y(x) = 
—b'(x). For large to only the first term on the r.h.s. of M .21 
contributes and one recovers essentially the bare disorder B. 
Decreasing to simply amounts to translate the curve upward 
along positive x, and yo increases as the curve x(y) cuts the 
axis x = closer to the minimum. It reaches it at the Larkin 
mass, solution of yo = Vc, i-e. 



TOfa 2 +^ lc = (ff7A) 1/(1+7) =--T c 



(7.33) 



For SR disorder yo = y c = e gives ml = g/e. Exactly as 



TO 



the solution acquires a cusp and one finds: 



b'(x) - b'(0) « y/-2(e-2()b'{0)x 



(7.34) 



i.e. the same result as (17. 19t . 

Although it is a priori not obvious how to follow this solu- 
tion for to < to c , the following remarkable property indicates 
how to proceed. If we compute the /3-function, i.e. the r.h.s. 
of (17.81 using 17.271 at to = m c and £ = £(7) we find that 
it exactly vanishes. Similarly the /3-function for b'(x) also ex- 
actly vanishes for all x > provided we use also J7.23l >. i.e 
all b'(Q) are defined as b'(0 + ). Thus at to = to c the func- 
tion has already reached its fixed point form x = x*(y), and 
freezes for m < m c . For the disorder correlators studied here, 
b(x) evolves according to 17.271 or J7.29I until m c where it 
reaches its fixed point b — b*(x), and does not evolve for 
to < to c . In particular y n — —b'(Q + ) freezes at m c and one 
has —md m yo = for to < m c , exactly as was discussed in 
the previous Section. 

The solution for m < m r is thus: 



x = (y/g) 1/7 + e 1 (y-yo)-T c m<m c 



-b'(0 + ) = yo= gf--y 



(7.35) 



where the parameter T c is defined in ( I7.33> . thus it exactly 
identifies with the zero temperature fixed point (17. 15l > with 
C = C(t)> as can t>e explicitly verified. This is easily un- 
derstood a posteriori, since the same functions appear and in 
both cases we have two conditions to fix the two undetermined 
amplitudes x(yo) — x'(yo) — 0. It does however heavily rely 
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on the exact power law form of the model, so it is not imme- 
diately obvious how it will extend to an arbitrary bare model 
B(z). One clearly cannot expect in the general case that con- 
vergence to the fixed point will be completed within a finite 
scale. The solution to this puzzle is given below. 

Similarly the solution for the Gaussian SR disorder correla- 
tor for m < m c is given by setting to = m c in ( I7.29l l. J7.30I 
with yo — e (which determines m e ). 

The result of this section thus provides unambiguously a 
solution beyond the Larkin scale which connects with the zero 
temperature fixed point. It justifies the previous Section and 
the value obtained for £. We found that for power law and 
gaussian models the freezing mechanism apparent in (I7.23> 
leads to: 



to > m c (7.36) 
to < to c . (7.37) 



- md m y(0+) = (e - 2C)y(0+) 
-md m y(0+) = 

The fixed point is reached at to = m c . 



D. Role of temperature 

In the case where disorder is relevant i.e. for 0(~f) > (i.e. 
2 < d < 4; d < 2 for 7 < j c = 2/(2 - d)) we found 
in the previous Section that temperature plays only a minor 
role since the convergence to the non analytic zero tempera- 
ture fixed point occurs on a finite (Larkin) RG scale. Whether 
it should be called a zero temperature fixed point can also be 
debated since it is reached when T m — T nic . A proper def- 
inition of the renormalized temperature may then include the 
denominator in (I6.2H . 

Let us now examine the marginal case #(7) = 0, 7 = J c (d) 
and d = 2 for SR disorder. We give here the main results, 
further details are examined in the Appendix|F| 

The analytic solution is given by ( 17.271 1 and yo given by 
f738l . where here f m = 4A d T/(e(2 - d)) does not flow as 
to is lowered. Let us examine the second derivative, 



fe"(0) 



~x'(yo) 



1 

72/o 



T 



-TO 



1 

e 

1+7 



(7.38) 



which is a rescaled version of ( 16. 1 31 . The first line in (17.381 
holds more generally (in the infinite UV cutoff limit) and to 
obtain the second we have set £ = C(t)> $ = an d as ~ 
sumed d < 2. Setting now 7 = j c , i.e. 8 = 0, we find that 
there is a transition at a temperature T = T C defined by 



<2-d)~ _ e(2-d), s 



4A d 



-T c = 



4A d 



(<h7< 



\ 1/(1+7) 



(7.39) 



such that for T > T c the solution is analytic for all to down to 
to = 0, given by 17.271 and b"(0) remains finite and given by 
(I7.38> . This is a line of analytic fixed points which terminates 
at T c . For T < T c the solution freezes as in the previous 
Section, and becomes non-analytic at and below the Larkin 



2 2C 



T r 



(7.40) 



The case d = 0, 7 = 1 corresponds to the logarithmically 
correlated disorder B(z) = —gln(a 2 + z). It has been studied 
for finite N in 1 30] where it was shown that there is a transition 
forany N atT c = yfg (g = a/N in the notations of Ref. 1 30]). 
The above result is in agreement with this value for T c . There, 
for N = 1, 2 there is also a line of fixed points for T > T c 
with a continuously varying dynamical exponent (and also one 
for T < T c with a different dynamical exponent and some 
form of RSB). Since the dynamical exponent is perturbatively 
related to b"(0), obtained above for infinite N, it would be 
particularly interesting to study the 1/N corrections in this 
case. 

Let us now examine the case of SR disorder ( I2.9I > in d = 2. 
More details are given in the Appendix |F| One has T m = 
2Th = (T/tt) ln(A/m). The analytic solution ( 17391 . ( TOOT l 

becomes 



x = hx(y /y) 
Va= gm 1 



e A -T/n 



(7.41) 
(7.42) 



with e = 2. Thus there is a transition at T = T c = 2ir. For 
T < T c , we find yo to increase as to decreases and reach yo = 
e at the Larkin mass. For to < m c the solution remains frozen 
to ( T7.41> with yo = e - F° r T > T c , we find that yo flows to 
zero and disorder is irrelevant. The physics is the same as the 
one contained in the variational method for the periodic model 
in d = 2 1 6] which exhibits a (so-called marginal) 1-step RSB 
solution. 

The case 7 > 7 c (d), (d < 2) is discussed in Appendix IEI 
Although an analytic solution exists as to — * and disorder 
is formally irrelevant, there are some freezing phenomena at 
intermediate to. It corresponds to the case where MP find, 
in addition to a RS solution, a 1-step RSB solution which is 
so called non-marginal (different in nature from the one step 
solutions obtained in the case 9 = 0). 



VIII. COMPARISON BETWEEN THE RSB 
AND THE FRG APPROACH 

In this Section we compare the FRG approach at large N with 
the GVM using RSB. Since the two methods study the same 
model in the same limit (large N) a precise connection should 
exist. 

We start by comparing the two methods at the level of the 
results of the calculations. We first perform the comparison 
for power law models. Then we generalize the FRG solution 
to arbitrary bare disorder correlator. Based on these results, 
we address the deeper connections between the two methods, 
and emphasize what we learn from them about the physical 
consequences. 
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A. Zero momentum correlation function from the 

FRG 

Our main result up to now is a non-trivial solution for the 
renormalized disorder correlator B{z) as a function of the 
scale parameter m, i.e. the effective action for the zero mo- 
mentum mode. Since this function is once differentiable, i.e. 
B(z) = B(0)+B'(0)z+O(z 3/2 ), we can extract from its first 
derivative the 2-point correlation function at zero momentum 
(see Section lirO : 



(v a (q) ■ "&(</)) 

G ab (q = 0) 



— (u a (q) ■u b (q / )) 
G ab (q)(2n) d 5 d (q + q') 



T 



= 0) = -^5 ab - 2 



B'(0) 



T A 
— — TrO ab — 



KM, 

2A d ' 



-d-2C 



(8.1) 



(8.2) 



where in the last equation we have used the definition J6.20I 
for the rescaled function b, and added the index m to recall its 
dependence on the mass. 

In the case 6 > 0, for the power-law models J2.10> . we thus 
find, using b' m (0) = — y from J7.28I 



GaMv = 0) = 2 5 (m 2C a 2 + f m J-~<m 
1 . i f e\ r & 
2A d y \-f 



d-2C 

-d-2Q 



(8.3) 



and, for m > m c the DR result 16.16k (16. 171 (where m c is 
determined by J7.33l ~) and, we recall, £ = e/(2(l + 7)). 

B. Explicit full RSB solution at large N 

Let us recall the RSB solution at large N and resolve carefully 
the MP saddle-point equation in presence of a mass. We only 
assume that there is indeed full RSB, to be checked a poste- 
riori. Let us first reexpress the general solution, valid for an 
arbitrary B, in a rather compact form. 

In the RSB method one first parameterizes the correla- 
tion matrix as G ab (k) = G(k, u) and the self-energy matrix 
TG~ b l {k) — (k 2 + m 2 )8 ab = a ab = <r(u), in terms of the 
overlap < u < 1 between (distinct) replicas a and b (and 
denote G = G aa ). The saddle point equations then read 



cr(u 

x(u c ; 

with 



= --B'(x(u)) 

= 2 f (G(k)-G(k, U )) 

Jk 



= 2T 



dw 



2Tcr / (w) 



k (k 2 +m 2 + [(t](w)) 2 



k 2 



(8.4) 



(8.5) 



(8.6) 



and E c = [er](u > u c ). The last two equations are the RSB- 
matrix inversion formulae; cr(u) is assumed to be continuous. 
Taking a derivative of J8.4> w.r.t. u gives 



a'(u) = c7'(u)45"(x(u)) 



1 



k (k 2 + m 2 + [cr](u)) 2 



(8.8) 



This equation admits two solutions: Either cr(u) is constant, 
or satisfies the marginality condition 



l = 4S"(x(u)) / 

Jk 



1 



(k 2 + m 2 + [er](u)) 2 



(8.9) 



We thus look for a solution of the full RSB equations (see 
Fig. |SJl with a non-trivial function cr(u) for u m < u < u c 
joined by two plateaus 



<t(u) = cr(u c ) , u > u c 

er(u) = a(u m ) , u < u, 



(8.10) 
(8.11) 



Similar forms are valid for G(k, u) and x(u). The breakpoint 
u c is related to the physics at the Larkin scale m c , which, at 
weak disorder, can be much smaller than the UV scale A, 
while u m depends on the IR cutoff rn. ( 18.91 also yields by 
continuity a closed equation which determines E c 



1 = 4B" 2T 



as well as 



1 



1 



k (k 2 + m 2 + E c ) J J k (k 2 + m 2 + E c ) 2 ' 

(8.12) 



1 = 4I 2 B"(x(u m )) 



(8.13) 



since [ct](u) = for u < u m . To solve these equations one 
firsts determines the function [er](u) (see below), then finds u c 
and u m . 

One can already note at this stage that j8.12l i is exactly the 
condition ( 16. 12t which determines the Larkin mass m c , equiv- 
alent to the vanishing of the replicon: 



i 2 — m for m < m c 



(8.14) 



and Y. c — (no RSB) for m > m c . 

To find [<t](u) for arbitrary B and cutoff, one notes 0] that 
with the help of d8.9l l and J8.4I > a can be expressed as a func- 
tion of [a] as 




(8.15) 

where (B") 1 is the inverse function of B" . Then one notes 
that u as a function of a of [a] is from d8.7l > simply 1/u = 
der/d[cr]. This yields immediately, using the chain rule: 



u = -4T 



k (fe 2 +m 2 + [o-](u)) 2 . 



fe (F+SP+I^MP 



M(u) = UCT(U) - [ dWCT(w) (8 - 7) xB '" [ {B " rl ( ^=w ), 



(8.16) 
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[o](u) + m 




FIG. 8: Full RSB solution for the function [er](u) + m and a fi- 
nite mass m. <r(u) has identical behavior, with two plateaus at val- 
ues <j(u = 0) and cr(u = u c ). In both cases, upon increasing the 
mass only, u m varies (increases) and the lower plateau moves up, the 
rest of the function being unchanged, see Eqs. 18.5 li f. in the text. 
The dashed line is the zero mass solution. The lower breakpoint u m 
reaches the upper one u c at the Larkin mass m = m c above which 
the solution becomes RS. The FRG gives exactly the lower plateau 
value for <r(u = 0) = a(u = u m ) and its m-dependence. From it, 
the full RSB solution can be reconstructed, see Section lVllIHI 



Upon inversion one obtains the exact function [er](u), and in- 
serting into ( 18 . 1 5i er(u). More precisely, we see that the sum 



is a m-independent function of u, and thus from 



fl8. 15i er(u) is also m-independent. Then one solves the self- 
consistent equation j8.12t for E c , and finally obtains u c from 
the above. The result can be written using ( 18 . 1 2i in the simple 
form 



-4T 



[Ik 



B'" (B") 



k w+^TF 



1 



- 4 Ik (k^+m'i)* ' ) 

(8.17) 

Thus u c depends only on the Larkin mass and is independent 
of m (See appendix|H]for another derivation and a discussion 
of this useful property). Similarly one obtains: 



I 3 

-AT-^-B'' 



(B'T 1 



Ah 



h = 



k (k 2 + m 2 ) 



2\Z 



(8.18) 



Let us apply these considerations to the power law model 
(12 . 1 Ot . For this model the Larkin mass is determined by 
( 17331 . Next one has: 



B'"{m-\y))=-l?^(^Y" . (8.19) 
4A d V 75 / 

In the limit of infinite UV cutoff A limit, using I 2 = 

(m 2 )~ e / 2 A d /6 and I 3 = (ra 2 )~ 1_e / 2 A d / '4 we obtain from 



M(u) = (Au) 
A 



1 



4TA d (l + 7 ) 
u m = m 9 1 A 
u c = m e JA 



(8.20) 
(8.21) 

(8.22) 
(8.23) 



with 6 = d - 2 + 2(, C = e/(2(l + 7)). Using (TQ31 one 
finds the m-independent result 



a(u) 



(8.24) 



In particular one has the value of the lower plateau (see Fig. 

ID 



cr m (0) = u(u m ) = 



2 - 6 u r 



-Am 



2-0 



(8.25) 



Let us already note the relation u rra 9 m 2(j m (0) = 1 which will 
be demonstrated to hold more generally below. 

C. Correlation function in MP solution compared 
to FRG 

The inversion formula yielding the diagonal correlation from 
the RSB solution is 



G aa (q = 0) 



T 



1 



<7 m (0) 



du H(u) 



u 2 m 2 + [cr](ii) 

(8.26) 

and is a sum of contributions from all overlaps < u < 1. In 
particular the contribution from states with zero overlap, i.e. 
the most distant states, is: 



G aa (q = 0)| u=0 =G{q,u = 0) = 



Ta m {0) 



(8.27) 



We can now compare with the FRG. One has, using 8 = d 

2 + 2C,2C-e(l + 7 ): 



G aa (q = 0)| u =o 



Ta m {0) 
m 4 

2(1 + 7) 



2 - 1 
1+27 



-T Am 



-1-6 



e 7 4A d (l+7) 

RG 



(7ff) 1+7 171 



= G F ^(q = 0) 



-d-2C 



(8.28) 



as given by (I8.3i . Thus, for this power-law model, we found 
that the FRG gives exactly and only the contribution from the 
most distant states (the lower plateau in the RSB solution). 
Before discussing the reasons and consequences, let us show 
that this feature is much more general than power law models, 
and holds in any case where full RSB holds. 

D. Solution of the FRG equation for arbitrary 
disorder correlator B 



In Section fVIICI we found how to continue the solution of the 
FRG equation beyond the Larkin scale. It involved freezing 
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of the to dependence of yo = — b'(0) at to = m c and worked 
only for two special forms of disorder correlators, which hap- 
pened to be already fixed point forms. It is important to find 
the solution for a more general form of the bare correlator 
B(z), and this is what we achieve here. 

Let us examine whether we can find a solution for any to of 
the FRG equation d7.8i in inverted variables 

md m x m (y) = {e-2()yx' m (y) + 2(x m (y)-y + yo , (8.29) 

which correspond to a more general function B{z). We take 
special care here to indicate that x m (y) is an m dependent 
function of y (we note x' m (y) = d y x m (y) and we recall that 
ym(x) = —b' m (x)). The idea is to play with the to depen- 
dence of yo = ?/o ( m ) since this is really all the freedom we 
have. Let us restrict our analysis for simplicity to T = 0, the 
generalization being straightforward. The definition of yo(m) 
is given implicitly by 



x m {yo(m)) = 
for all to. The total derivative thus vanishes: 



(8.30) 



m—(x m (y (m))) = md m x m (y ) + x' m {y )md m y 
am 



= 



(8.31) 



Together with ( I8.29> at y — yo, it yields (recall that x m (yo) = 

0): 

(md m y + (e - 2C)yoK„M = . (8.32) 
There are only two possible solutions: 



md m y + (e - 2C)y = 
x'{y ) = 



(8.33) 
(8.34) 



The first holds before the Larkin scale and the second, which 
implies a non-analytic b(y), beyond. We now want to find the 
solution beyond the Larkin scale, i.e. assuming that x' m (yo) — 
0, together with x(yo(m)) = 0, which of course implies 
md m x m (yo ) — 0. 

Equation ( I8.29> with yo = yo (to) is trivially separable and 
admits the general solution 



x m (y) = m 2C § ^ 



V 



,2C 



1 

00 Am' 



/ ^ yo ( m ')m'- 2 < , (8.35) 



where for now yo(m) is arbitrary and so is the function 
4>(j/). (It will be identified below with (— B')^ 1 as in Section 
IVIH . The first condition one must impose is the definition 

x m (yo{m)) = 0, i.e. 



(8.36) 



AA d m 2 ^ t J e' 
/ — rVaim )m 



which should be valid both for m > m c and m < m c . Taking 
TO dm °^ ^-36> yields, using J8.36I I again 

I i JUL®' ( ^o( TO ) 
e + 4A d UA d m^-* 



[md m y + (e - 2()yo] = . 



(8.37) 

In order to satisfy this equation, at least one of the factors 
must vanish. The regime m < m c corresponds to the first, the 
regime m > m c to the second factor being zero. 

For to > m c one has \md m yo{m) + (e — 2()yo(m)] = 
leading to 



yo (m) = Am 2 ^ £ 
and the above solution becomes: 

y \ l 



= TO 2f $ 



AAdm^-t 



(y - 2/0 ) 



(8.38) 



(8.39) 



This can clearly be identified with the analytic solution of 
the self-consistent equation i7.1\ found before in Section lvTTl 
and thus implies that $ is the reciprocal function of —B'. 
Eq. J8.36l > is trivially satisfied by 



\AA d m 2 <- 
Applying —B' to 18.401 fixes A to be 

A = -AA d B'{0) 



. 



(8.40) 



(8.41) 



and one recovers the dimensional reduction result. 

The interesting new information is obtained for m < m c . 
Then the first factor in d8.37t vanishes, i.e. 



e ^ AA d \lA d m*-* 



(8.42) 



Deriving (I8.35i w.r.t. y one sees that J8.42I correctly implies 
x'{yo) = , (8.43) 

thus the solution for b'(x) has a cusp. J8.42l > determines the 
function yo(m) for m < m c . Note that if the power law in 
the correlator holds only asymptotically, yo(m) will nicely 
converge to a constant (for the right choice of Q due to the 
asymptotic power law tail, but may vary arbitrarily according 
to the irrelevant corrections to power law. This is studied in 
more details in appendixlHl 

It is convenient to rewrite the final result, i.e. Eqs. J8.36i . 
(I8.42> in the form: 



b' m (0) = AA d m 2 ^B'(x m {0)) 



B"(Xm(0)) 



(8.44) 
(8.45) 



^C(0)to'- 2 <, (8.46) 

TO 



where we use the notation X m (Q) 
nection with the RSB solution becomes obvious in this form 



X™ G (0). The con- 
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Comparing with J8.13> . the equation J8.45> of the FRG solu- 
tion identifies with the marginality condition at u = u m , the 
lower plateau of the RSB solution, see Fig. [8] It allows to de- 
termine Xm(0); the two other equations are self-consistently 
obeyed and give b' m (0). Comparing with J8.4i at u = u m 
yields the identification 



Note that Eqs. ( 18. 511 and J8.52i have been hypothesized by 
Parisi and Toulouse for the SK-model |61]. However, it has 
been shown that there they are only approximately satisfied, 
see e.g. ill H. 

Analysis of these equations shows that, up to the break- 
point, one has: 



Xm(u = 



-- 0) = x(u = u„ 
(0) Ta(u m ) 



X™ G (0) 
-^(0). 



-d+2C 



2A H 



(8.47) 
(8.48) 



and thus we obtain: 

J («z)=G GVM feu 



a m {u) 


= CT (w) 


U> U m 


(8.53) 


<J m ] (u) + m 2 


= [<ro](u) 


u > u m 


(8.54) 


<J m {u) 


= fTm(O) 


u < u m 


(8.55) 


[Um] (u) 


= 


u < u m . 


(8.56) 



^.FRG i 



0) — qGVM (q\\_ Q (8.49) u m is thus uniquely defined from the solution at zero mass by 



It thus holds for an arbitrary disorder correlator, provided a 
solution to Eqs. J8.44I . J8.45> exists, i.e. for the class of func- 
tions B(u) which yield full RSB (also called continuous RSB) 
within the MP approach. Of course, equations J8.44i . J8.45i 
were derived without any assumption about replica symmetry 
breaking. 

Extension to T > is obvious. Adding the last term of 
A7.8I and following the same steps as above, one finds: 



md n 



yo 



2C)yo 



(8.50) 

Vanishing of the first factor yields the finite T analytic solu- 
tion studied in the previous Section (equivalent to the RS so- 
lution of MP). Continuation beyond the Larkin mass implies 
x' m (y o) — 0, in which case the additional temperature term in 
(17. 8> vanishes and one is back to the T = equations ( 18.451 1. 
(I8.44> : Thus only the value of the Larkin mass depends on 
temperature, everything else is independent of T. 

E. Full RSB solution from the FRG result 

In the previous Section we have shown that the FRG yields 
<7 m (u = 0) (via b' m (0)) i.e. the value of the RSB function 
only at u = 0. In fact, as we now discuss, by varying the 
mass one can scan the whole function <r(u) of MP for any u, 
and thus the FRG yields the same information as contained in 
the function <r(u). Remarkably, we can obtain an explicit ex- 
pression for cr(u), even though the argument of this function, 
the "overlap" is not obviously related to any quantity in the 
FRG. Furthermore, we can also compute the full correlation- 
function of Mezard Parisi, if one knows only <j m (0) for all m, 
which is given by the FRG. 

Thus from now on we assume that we know only a m (u = 
0) as a function of m through the FRG, together with some 
general properties of the MP solution. As we have already 
found in Section IVIIIBI and is shown more directly in Ap- 
pendix [G] the GVM saddle point equations, upon assuming 
full continuous RSB, satisfy the two "RG equations" 



d m a m (u) = 
d m {[a m ](u) +m 2 ) = 



(8.51) 
(8.52) 



valid for any u such that cr'(u) 7^ 0. One can thus relate the 
solution <7 m (u) at finite m to the solution at zero mass ao(u). 



(T m (0) = CT (u m ) 
[era] (u m ) = m 2 . 



(8.57) 
(8.58) 



Indeed one has, taking derivatives of d8.57i and ( 18.581 1 w.r.t. 

2 

m . 

<9 m 2(7 m (0) = a' (u m )d m 2u m (8.59) 
u m o-o(u m )9 m 2u m = 1 , (8.60) 

where here we introduce for convenience d m 2 = -r^d m . 



These two equations give 



1 



(8.61) 



d m 2a m (0) 

One thus finds that the function oo(u) is implicitly given by 

1 



cr U 



<9 m 2(7 m (0) 



MO) • 



(8.62) 



Since er m (0) can be extracted from the FRG, we see that we 
can obtain the full function a (u) from the FRG. 

One notes that the upper breakpoint = u c is indepen- 
dent of m. As shown in Section dVIIIBL u m increases upon 
increase of m, and reaches u c at the Larkin mass, i.e. for 

lim e = U c . 

Let us show how one can recast the correlation function of 
MP, given in ( 18.261 at zero momentum, entirely using FRG 
data. 



G aa {q = 0) 



1 



fm(0) 



rn 

du 
— 2" 



u 



1 



[MM 



Using our previous results gives: 



G aa (q = 0) 



1 



<Tm(0) 



m 
du 



1 

Uc. 



(8.63) 



(8.64) 
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M(u) 

ti 2 im 2 , 
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Shifting from the variable u to the variable fi defined by 
u M = u, one finds that the correlation function can be ex- 
pressed entirely from the knowledge of er m (0). To see this, 
note that 



du 



-d 



= -d(S m 2a m (0)) = -d 2 m2 a m (0)d(m 2 ) 



(8.65) 



This gives 



G aa (g = 0) 1 | <r m (0) 

i 2 to 4 



T 



w, r /V 771 mi 



(8.66) 



where we have used that [ao](i 



/i 2 . After an integra- 



tion by part and using again J8.6H at u = u c one finds the 
remarkably simple formula 



G aa (q = 0) 
T 



dw 2 „ , , 1 
Ji-0 ^(0 + -3 

/J, 4 777^ 

da M (0) 



1 



(8.67) 



valid for m < m c . Recalling the relation between cr m (0) and 
-B'(O) we see that the MP result is a simple average of the cor- 
relations corresponding to masses between m and m c . Note 
that one can derive a similar formula for G(k, u) obtained by 
MP as a function of cr^(0). 

In summary, although strictly speaking our FRG result 
gives only the contribution of distant states to the 2-point cor- 
relation function, it does allow to obtain the whole MP result, 
although we do not yet have a derivation within the frame- 
work of the FRG alone. One should also note that the formula 
(I8.46> is in a sense equivalent to the inversion formula of hi- 
erarchical matrices which relates x(u) = 2 J, G — G(k, u) to 
the self energy cr(u). This raises the question of whether the 
FRG equations "know" about ultra-metric matrix inversion. 
These result hold for continuous RSB and the case of non- 
marginal RSB (when the marginality condition is not obeyed) 
is discussed in AppendixlEl 

F. Discussion: Explicit versus spontaneous 
replica symmetry breaking 

Let us examine what has been achieved and how it compares 
with other works. 

We are interested in the behavior of the effective action of 
the replicated field theory for N large. Let us focus here on 
the uniform configuration, for which T(u) = L d T{u), where 
we denote T(u) the effective action per unit volume. TY(u) 
then represents the free energy per unit volume, depending 
on a set u a , a — 1, . . . , n in presence of external sources J a 
which impose field averages (u a {x)) — u a . The usual free 
energy T is recovered for u — 0, T — TT(u = 0). 



The saddle point method allows to write in the limit of infi- 
nite N: 



e -i» = C -L d r(u) 



E 



(8.68) 



where f (v) satisfies the saddle point equations j4.2H4.3t . and 
J4.5i . and the denotes a sum over saddle points whenever 
more than one solution exists |75|. To be accurate the saddle- 
point method computes 



lim : 

N^oo NL d 



T(u = vVN) 



(8.69) 



where the limit is taken at fixed v. 

As mentioned in Section HV Al the saddle point equations 
J4.5i contain both the FRG and the GVM. They depend on 
the set of v a , and when expanded in cumulants, taking all 
v a i, = v a — Vb 7^ 0, they lead to the FRG equations. This 
approach clearly consists in imposing an explicit breaking of 
replica symmetry. Also we expect that in that case a single 
saddle point exists. This is indicated by the fact that the quan- 
tity B'^v 2 )^ 1 = S"(xu) _1 — 4/ 2 plays the role of a replicon 
eigenvalue and remains frozen and positive for v > 0. 

By contrast, MP studied the case where all v a b = and 
found spontaneous replica symmetry breaking, i.e. multiple 
saddle points, differing by permutations it of replicas. 

We can now make contact between the two approaches and 
understand why we have obtained via the FRG the correla- 
tions of MP corresponding to the distant states. Let us focus 
on the mode q = 0, and define the center of mass variable 
u := yjii(x) (i.e. without rescaling in N) and consider: 



-L d Y: a JaU a 



Z(J) = / dui . . . du n Pv(ui) . . . Pv(u n ) e 

(8.70) 

Py('u) is the probability distribution of it in a given disorder 
configuration: 



P v (u) = / V[u}5 u 



1 

1? 



u(x) 



-Hv[u]/T 



(8.71) 



In the present paper we have computed J8.7QI >. scaling J a ~ 
V~N, and taking all J a different in order to impose all v a b ^ 
~ O(l). Because of this scaling with N we obtained a 
different saddle point than MP (shifted by j, see ( I3.12» . and 
since all J a are different, this saddle point has explicit RSB. 
According to Section||ll]this gives T[u] when u scales as yN, 
i.e. we determined the averaged probability ( 12.341 1. ( 12.351 1 for 
fixed w = u/ \/N. 

On the other hand MP found that: 



P v {u 1 )...P v {u n ) = Y, e-^ G ^^= >^ 



(8.72) 



One should in principle be able to recover the 2-point corre- 
lation function d8.26> obtained by MP adding small sources 
J a as in (I8.70> and taking derivatives at J a = 0. Clearly, 
to reproduce the MP result, these should be taken as J a = 
O(N ) — > and not 0(vN) so as to maintain the unper- 
turbed MP saddle point. For instance the diagonal 2-point cor- 
relation function is obtained using J\ = J, all others J a = 0, 
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and differentiating twice w.r.t J. (The off diagonal one in- 
volves ( J, — J, 0, . . . )). Equivalently it should be obtainable 
from the effective action for u a t, = u a — Ub = 0(N°). Thus 
scaling u a b ~ V~N and J ~ >/N as was done here selects 
the distant states in the MP solution. The fact that we obtain 
exactly the MP result for these states shows that there is no 
intermediate scaling regime. 

We emphasize that our primary aim here is not to recover 
the MP result, but to understand what exactly the FRG pre- 
dicts, in view of getting a better understanding of FRG within 
e.g. the e-expansion. Extension of the FRG beyond the Larkin 
scale requires giving a meaning to the limit u — > + . We find 
here that what the FRG actually computes (from b'(0 + )) is 
a second moment of w in presence of a small extra external 
source vNj a such that all v a b ^ 0, i.e. an average, such that 
when there are several states, the different replicas are chosen 
in maximally separated states (u = 0). 

Note that the quantity computed by the FRG specifies the 
system's preparation, while such a procedure still has to be 
worked out for the MP solution. In presence of a broken sym- 
metry this is an important issue, and the FRG gives a natural 
solution. 

It would be interesting to understand further the limit 
Vab —> coming from our solution, which one can call the 
crossover from FRG to RSB. It is clearly non-trivial. For 
instance, one question is what we get if we take a source 
J = (jy/N, — jV~N, 0, . . . , 0) so that we still have sponta- 
neous RSB in n — 2 copies, or if we divide the replicas in two 
groups of n/2 each, J = (jVW, ..jVN,-jVN,0, ..-jVN) 
so that RSB persists within each packet. 

Another important issue is what happens at large but finite 
N. For any N, if one parameterizes the 2-replica part of the 
effective action using R(u a b) = NBpf(u^ b /N), one can write 
the 2-point correlation function as 

G ab {q = 0) = -2B' N {0)m- 4 (8.73) 

for a b. We have determined the function Bn(x) for 
x = 0{N°), i.e. u ab = 0(y/N). To obtain B' N (0) for 
finite N however, one needs a priori to know R(u a b) for 
u ab ~ O{N ), i.e. B N (x) for x = 0{1/N). The two could 
be the same, or there could be a boundary layer of size 1 /N, A 
priori the knowledge of this requires including the 1/N cor- 
rections in the FRG equation (as is examined in 1 47]). This 
may help to better understand the connection of this regime to 
RSB. This is important since there are cases (e.g. for d = 0, 
6 < 0) where we know that Paris-type RSB cannot survive at 
finite N. 

G. Interpretation: Comparison with BBM approach 

A previous study |42] aimed at connecting the RSB solution 
to the FRG. The authors defined, for each configuration of 
the disorder, an "effective random potential" U~v(4>o) for a 
given mode (e.g. the zero mode). Starting from the MP so- 
lution ( 18.721 1. they computed the second cumulant of Uy and 
showed that it exhibits a non-analyticity, reminiscent of the 
one found in the FRG. A parallel was drawn with a d = 0, 



N = 1 toy model where Uy satisfies an exact RG equation of 
the Burgers-KPZ type with random initial conditions. Shocks, 
well known to develop in this equation, provide an appeal- 
ing physical picture for the singularities in the energy land- 
scape responsible for the non-analyticity in the FRG beyond 
the Larkin length £lo ■ 

Comparison between this study and the present one shows 
several important differences, with interesting physical con- 
sequences. The scaling in L, in N, and the definition of the 
"renormalized" disorder are different. As here, the authors of 
Ref. l42Tl focus on the zero mode, but with a different scaling 
with system size: They define 0o : = L~ d l 2 J u{x) = L d / 2 u, 
such that fluctuations of 4>q remain of order one. Other quan- 
tities are: 

U v (<h) = -TbxP v (<po) (8-74) 
Uv{4>o)Uv{4>' ) = -Rbbm(0o - 4>'q) = B B bm((4>o - 4>'o) 2 ) ■ 

(8.75) 

This definition means that the mode <f>$ sees an equivalent 
d = toy model, with random potential Uv{4>a)- Compar- 
ing with d2.35i and ( I2.37K we see that since the rescaling in 
L is different there can be no relation between -Rbbm (0) or 
B BBM (0) and the two point correlation function, neither the 
one of MP, nor the one obtained here in the FRG. To obtain 
the 2-point correlation of MP one would still have to solve the 
toy model defined by Uv(<f>a), i.e. compute: 

J d0j . . . r <t>l ■ (8.76) 

This task is difficult, since it requires not only the second cu- 
mulant, but also higher ones (not computed in 14211 '). More 
importantly, it requires the large argument behavior c/)q — 4>' of 
i?BBM (y)> n °t obtained in 1 42], were attention was focused on 
the Larkin scale (see below). Thus the information contained 
in i?BBM is physically interesting but not obviously related to 
large scale correlations. It is in a sense (e.g. for the d = 
case discussed in I42l0 closer in spirit to Wilson-Polchinski 
type RG l62ll versus an RG based on the effective action (see 
145] lH). 

Reexpressed in the variables of the present work, the result 
of 14211 reads: 

/ L \ d/2 

B{u 2 ) = B'{0)u 2 + c — u 3 , (8.77) 

where c is a constant. Because of the different rescaling, the 
non-analytic term has a coefficient growing with the system 
size, which expresses again that it is not an effective action. 
However, since the u 2 term is simply the bare disorder, and 
the non-analytic term involves only the Larkin scale £lo> it 
seems that this carries information for and only for the physics 
below and around the Larkin length, and does not contain any 
information about large scale behavior. Thus, despite exhibit- 
ing shock behavior at the Larkin scale, we think it has little to 
do with the FRG as a perturbative method to obtain large scale 
behavior. Not surprisingly then, -Rbbm is non-perturbative in 
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e = 4 — d, contrarily to the one obtained in standard FRG, 
which is of order e. 

Another important difference with the present approach is 
the scaling with N. The approach of [42] used the unper- 
turbed MP saddle point and thus, as was extensively discussed 
in the previous Section, it focuses on u a i, = O(N ) while we 
focus on u a b = O(N), (i.e. v 2 ab ~ 1/N there and v^ b ~ 1 
here) . Further work is needed to connect these regimes. On 
the other hand it seems that the thermal boundary layer can be 
found within this approach |77||. 

IX. DISCUSSION AND CONCLUSION 

In this paper we have studied the FRG at large N. From an ex- 
act saddle point calculation of the replicated effective action 
at large N we have derived the exact renormalization group 
equation, valid in any dimension d for infinite N, for the field 
theory of pinning. It is expressed as the /3-function for the 
second cumulant of the disorder correlator, and is exact as the 
second cumulant satisfies a remarkably simple closed equa- 
tion. To order 0(e) it agrees with the one derived by Balents 
and Fisher. 

This result teaches us a lot about how the FRG works and 
helps put the FRG approach to the e expansion on more solid 
grounds. Since here the FRG flow equation is formally equiv- 
alent to a self-consistent saddle point equation, it is fully inte- 
grable, i.e. one can follow the RG flow from any initial con- 
dition. It is thus possible to examine in detail what happens 
around the Larkin length and how the disorder correlator de- 
velops the non-analyticity. Let us emphasize that this is the 
first time that the emergence of non-analytic behaviour in the 
FRG is proven rigourously, beyond perturbative calculations. 
Indeed, the 1-loop FRG is insufficient per se to provide such a 
proof since the runaway of R""(0) could very well be argued 
to be the analog of the famous Landau ghost, i.e. a flow to a 
strong coupling fixed point without the need, or better the pos- 
sibility, for renormalization within a non-analytic functional 
space. Here we demonstrate that this is not the case, at least 
for infinite N. 

If we had restricted the analysis to the self-consistent equa- 
tion, the continuation beyond the Larkin scale would have 
seemed quite problematic. Remarkably, the FRG equation, 
equivalent below the Larkin scale, provides an unambiguous 
way to continue the flow equation beyond the Larkin scale. 
Even more remarkably, its solution reproduces exactly the 
small overlap result of the full RSB solution of MP, a non- 
trivial result which, within MP, cannot be obtained without 
constructing the full RSB solution. The mechanism for this 
seems to be that the FRG solution in that case naturally satis- 
fies the so-called marginality condition. In fact, it turns out to 
be equivalent to it, and we were able to find a formula yield- 
ing the complete RSB solution for all overlaps. This is striking 
since we did not make any assumption about Parisi RSB. We 
avoided the issue altogether by using a method where RSB is 
not spontaneous, but explicit. 

Given that the validity of the Parisi Ansatz, e.g. for the SK 
model, has not yet been proven (despite recent progress ll63lo . 
it is interesting to find a method which does not rely on RSB. 



In fact there may well be deeper connections to be unveiled 
between the Parisi algebra of ultrametric matrices and the type 
of singular differential equations arising in the FRG. Another 
example where a RSB solution can equivalently be obtained 
via an RG type equation is the Derrida Spohn solution of the 
DP on the Cayley tree l3lll . This has inspired a similar solu- 
tion for a model with translationally invariant disorder in 1 30]. 

We have thus shown agreement with the main results of 
the full and the marginal one-step RSB solutions of Mezard- 
Parisi. This is also interesting since it has been widely debated 
J64ll whether the RSB method captures the physics: Our re- 
sults raise no doubt for infinite N. 

More puzzling is the situation for SR disorder. There MP 
find both a stable replica symmetric solution and a 1-step so- 
lution where minimization over the breakpoint has to be en- 
forced (marginality condition violated). For gaussian disor- 
der both solutions of MP have ( = (2 — d)/2. Similarly the 
FRG naturally finds the finite temperature RS solution with 
£ = (2 — d)/2 (and one fixed point solution at T = with 
C = 0). A non-analytic solution also seems to exist in the 
FRG, and work is in progress to analyze it further and eluci- 
date whether it is related to the non-marginal 1-step solution 
of MP. Let us note that in d = we essentially know that 
(apart from 7 = 1) RSB does not hold at finite N (the phase 
transition predicted by the GVM in that case must go away 
at finite N,T > 0). So there is little doubt that the correct 
branch at finite N is the RS one, as also given by the FRG. 
For the DP problem with d = 1 on the other hand, it is not yet 
clear whether both branches (a T = fixed point starting from 
( = or a finite T non-analytic solution with £ = (2 — d)/2) 
can coexist. One scenario is that they would cross over at 
some lower value of N — N c yielding the upper critical di- 
mension of KPZ. The calculation of the FRG /3-function to 
next order in 1/N should shed light on this question, and is 
thus of high interest. It is presented in l47ll . Our method thus 
provides a unique candidate for a field theory of the strong- 
coupling phase. 

To summarize, the present method is promising in solving 
mean field models, by using explicit rather than spontaneous 
RSB. It would be of interest to investigate whether other mod- 
els like the SK-model could be solved via the same route. 
More importantly, it may open an alternative road to tackle 
disordered systems from a different direction than expanding 
around RSB saddle points, a task which still has to be accom- 
plished. Of course, in the end, the same difficulties may be 
in store. They could hide in the subtelties due to the non- 
analytic behaviour of the /3-function at large N. However we 
are optimistic, since we have understood the infinite- A limit, 
at least in the full-RSB case. Also, a solution has been found 
for N = 1 to two II43L l44l I.55I1 and three loop order |65B, and 
for finite N at 2-loop order. 

Let us close by indicating that many extensions of this work 
are possible and some in progress. One example is the random 
field problem, still under intense debate 1 66 , 67 ] , for which we 
have also computed |68] the effective action at large N, and 
at 2-loop order. Finally, the same method applies to the dy- 
namics, classical or quantum: it has been shown |9] that the 
mode-coupling approximation in glasses |69] identifies with 
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(non-marginal) mean field (large N) dynamics, exhibiting ag- 
ing solutions. However this picture leaves out thermally ac- 
tivated processes, and our 1/N method may be promnising 
there too. 
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APPENDIX A: VARIATIONAL 
FORMULATION 



Let us extend the variational method of Ref. 1 24] to the case 
where the average of the field u a {x) is fixed to a non-zero 
value denoted here u a (x) — \/~Nv a (x). One defines the vari- 
ational Hamiltonian and free energy: 



H 



\arl u ] — 2 



[U a {x) - U a (x)} 



3'!/ 



x(Gvar)za,y& [u a {v) ~ U b (y)] 
^"varlGvar, V] = -Tin / V[u] exp (-7Y V ar M/T) 



(H-H 



var / ■}-[ 



(A.l) 



which satisfies (for n positive integer) the usual bound T = 
-TlnZ < .F var . Here H = NS[u,0] defin ed in t he text, is 
the replicated Hamiltonian. Comparing with ( I3.27> one finds 
that 

var [G var ,v]/T = T [G var ,V,U], (A.2) 

where the last argument indicates that for TV finite U(x) 
should be replaced by U(x)\ in the infinite- N limit U = U. 
Restricting to a bare model with only a second cumulant one 
finds (omitting the bars on v): 

U(vv(x), (G v ) xx ) = J2 B(v ab (x) 2 , {G v )t) 

ab 



B(v 2 ,z) = 

Q aa 



d N w 



(2n) N / 2 
G^ 3 G a ^ G 



w2 ' 2 B(v 2 



n v-w _ ur 

(A3) 



In general, B(s, z) is a function of two variables, which be- 
comes a function of the sum B(s,z) — > B(s + z) only as 
N — > oo, since in that limit (vw) = and (w 2 ) = N, without 
fluctuations, in the gaussian measure ~ c~ w I 2 . 

APPENDIX B: EFFECTIVE ACTION IN 
NON-UNIFORM BACKGROUND: GENERAL 
FORMULATION 

In some applications bilocal terms may already be present in 
the starting action. Let us thus give a more general and com- 
pact result, which also includes that case. It is derived by a 
simple extension of the methods of SectionllTll 



Let us consider a TV-component field <p x , with components 
4> % x , i = 1, . . . , N, which can carry other indices, coordinates, 

or be a set of fields, etc A general 0(7V)-invariant form 

for the action functional is 



S[cj>] 



N 



(B.l) 
(B.2) 



where S- m t [ip] is a functional of the bilocal field ip xy (which is 
also a bi-index matrix if the field <j> carries other indices, etc.). 
Then the effective action associated to S can be written as 

m = \(b:C- 1 :cb + NT ^]+T 1 [ij} + ... , (B.3) 

where T°[ip] is a functional of the bilocal field ip xy in ( IB.2i 
and satisfies the self-consistent equation: 



ST° 6S int [ X ] 

m - 



Slpxy^' SXxy 

Xxy = Ipxy + G[lp] X y 



G[ip] xy = [C- 1 + 2 



Xxy is another (set of) bilocal fields. 



Sip 



(B.4) 
(B.5) 
(B.6) 



APPENDIX C: CALCULATION OF HIGHER 
CUMULANTS 

In this Appendix we compute the third and fourth renormal- 
ized cumulants of the disorder. One uses the parameterization: 
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The equality of ( 14. 1 31 (pushed to the fourth cumulant, i.e. the 
above formula) and of (14. 14i using d4.20i . implies: 
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Thus we need compute the terms with one and two free replica matrix M, but only terms with zero, one or two replica sums, 
sums, Xab and Xab- Because of it means that we need The expression of (M 2 ) ab given in <4.25> is thus sufficient, 



to compute Xab' Xa\ Xg h_ aR d X<?' '■ To compute them we 
use the definitions (I4.12> . ( I4.18> . We thus need powers of the 
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and we also need: 
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dropping all terms with three or more sums. One then finds: 
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which yields Xab anc * Xab usm g < M-.24i . 



1. Third cumulant 

To obtain the third cumulant we now insert jQh m EST One can rewrite B"(x a b)' indeed taking the derivative of ( I4.29i with 
respect to v 2 b shows that: 



Kb 



{Xab) l + ihKb 

This becomes, regrouping the terms in S' labg and dividing by the common denominator 1/(1 + AI2 B" b ) we obtain: 

^^S'labg + ^l2 B ab(yS'l,aag + 7^^'l.bbg) 
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This first yields: 



4TI 2 B"(0) 
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1 + 47 2 S"(0) 
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Inserting this back yields: 
S" 



IT I2B 

l-"'-c - T~~md77^( B 'ac + B 'bc)+ 8I 3 B ab 



l + 4/ 2 B"(0) 



(s; 6 - s'(o))(s; c + s^j - b; c j% c + l -{B' ac ) 2 + l -(& bc ) 
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FIG. 9: Graphical representation of the fourth cumulant. The notation is explained in |47]. Each diagram corresponds to a square bracket in 
Eq. <C2U . in the same order. 



In terms of functions it gives: 



S[(x,y,z) = 



2TU 



l + 4/ 2 S"(0) 
+%hB"(x) 



B"(x)(B'(y)+B'(z)) 



(B'(x) B'(0))(B'(y) + &{z)) B'(y)B'(z) + \{B' {y)f + \{& (z)f 



(C.15) 



Integrating once, this yields the simple expression for the third cumulant given in the text. Note that, up to terms which vanish 
at n = 0, it can be expressed only in terms of the function B'(x) — B'(0). 



2. Fourth cumulant 



From iC.5\ one has: 



^ S'l,abe — 22 — 

e e 

^2 Q'labef = X] [t 

ef ef 



B ab 



where we have used B"'(xf l j) — B'" b j ((1 + 4l2B" b ) 3 obtained by further differentiation of ( 14.291 with respect to v\ b . We also 



Kb 



-Y, 



abe 



-z, 



B 



abef 



ab 



ab 



2(1+4/^)3 



^abe^abf 



(C.16) 
(C17) 



define 



Yabe — TXab 



and 



Zabef = T Xab 



4/ 2 
= 2h 



2^haae + ^l,bbe) + ^l.abe 
- -j^B'aa + B 'bb)( B ae + B' be ) 



2I 2 T(B' ae + B' b 



be) 



^{B'aeB'ae + B 'be B be) 



B 'ab( B ae + B 'b 



be I 



-81, 



+ 8/4 



B 'ae B be 



B be B bf) 



(C.18) 



abef ~ 2^1,aaef ~ ^M&e/J ^ 2 ^ X^l.ae/ + ^l,bef) + ^hT(B' ae B' af 
~ ( B 'aaS'l,aef + B 'bb^'l,bef) ~ ( B 'ae^'l,aaf + B 'be^'l,bbf) + ( B 'afS'lJae + B 'bfS'l,fbe) 
+ B 'ab(S'l,bef + ^l.aef) + ( B 'ae^'l,abf + B 'be^'l,abf) ~ ( B 'af^'ljbe + B bf^X,fae) 
B 'ae B 'ae( 2B 'af + B 'ef) + B be B be( 2B bf + B 'ef) ~ ( B ae B af + B 'be B 'bf) B 'ef " 3B 'aa( B 'af B ae + B 'bf B 'be) 
~ 2B 'ae B 'be( B af + B 'bf + B 'ef) + 2B 'ae B 'bf B 'ef + 2B 'ab( B 'af B ae + B 'bf B 'be + B 'ae B bf) > 



which should be further symmetrized over e, / for later use (indicated by sym e j below). 

The fourth cumulant equation yields, regrouping and simplifying the denominators (and using also the third cumulant equa- 
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tion): 
1 



'labef 



B be B bf) 



ab ab 
{ B aa^l,aef 



■2I 2 T S ym ef [B^ b (S' ltaef + S' 1<bef )] 

B'bbS'l,bef^ ~~ ( B ae^l,aaf + B be^X,bbf) 



+4I 3 TB': b (B' ae B' af 

+ ( B 'afS'lJ'ae + B 'bf^'ljbe) + B 'ab(^'l,bef + S' lae j) + (B' ae S' lab f + B' be S' lab f) — (B' a fS[j be + B'bfSljae) 

+8/ 4 S" b sym e/ B' ae B' ae (2B' af + B' ef ) + B be B' be (2B' b} + B' ef ) - (B' ae B' af + B' be B' bf )B' ef - "5B' aa {B af B a 

- 2B ae B 'be( B 'af + B 'bf + B ef) + 2B ae B 'bf B 'ef + 2B 'ab( B af B ae 

Setting a = b and solving one finds: 

8T 2 / 2 5"'(0) 



B bf B L + B 'ae B 'bf) 



1 ~, _ 4I 2 TB"(0) ~, 
2 Qaaef ~ l + 4/ 2 B"(0) Syme/ W 

This gives the final result for the fourth cumulant: 



SI 3 B"(0)T 
(l + 4/ 2 5"(0)) 3 T 1 + 4J 2 B"(0) 



B 'ae B af 



B 'bf B 'be) 
(C19) 



(C.20) 



^abcd 

48 
+192 



= Sym, 



abed 



-^ B B ab B ac B ad + ^( B ab) 2 B ac B ad + 2 ( B ab) 2 B 'ad B 'bc - iB ab B ac B ad B bc + B 'ab B 'ad B 'bc B 'cd 



h 



4 (- B q) B 'ab B 'cd B ac ~ iB 'o( B 'ab) B 'cd B ac 



& B B ab B ac B cd B ac 



4( B 'ab) 2 B 'ac B 'cd B ac 



" ^ B 'ab( B 'ac) 2 B cd B a, 



Q zV JD f £?' £?' jd" q tV /V ryt r>f r>lf q zV r>/ r>f ryt fytt i a ( ryt \2 ryt ryt ryft , a ( ryt \2 ryt ryt ryft 
-^ B B ab B bc B cd B ac ~ ° B B 'ac B bc B ' cd B ac - & B ab B ac B bc B cd B ac + ^\ B ac) B bc B cd B ac + A \ B ab) B ac B ad B ad 

- 4:(B ab ) 2 B'^B'^B 1 ^ + 2B'„kB'„„B'kjB'„jB"„a + ( BLiS 2 ( B'„,i) 2 B 1 !. 

-( B ad) 2 ( B cd) 2B bd 



J 1J ac 1J bc 1J cd 1J ac 

o r' R' R' r' r" 
2B ab B ac B bd B cd B ad 



J 'ac be 
( B ab) 2 ( B cd 



A{B' ) 2 B' ac B' cd B bc - AB' (B' ac ) 2 B' cd B bc 



+192T 

+32T 
+48T 2 
-128T 2 



O T2 1 T2 1 £?' f?" 
ZB ab B ac B ad B ab 

+ ( B abf ' B cd B ac 



B 'ab B 'ac B 'ad 



ab n ad n bc n ab - An n ab B ad B ac 



% B 'ah B 'ad B 'hr B "h ~ ^ B Q B 

hh 



( ryt \2 ryt ryf o £?' £?' T2 1 fi" 9 £?' £?' 
l^W B ad B ac ~ 2B ab B ad B bc B ac ~ 2B B ab B cd B ac 



ZB ab B ac B ad B 



1 + 4J 2 B&' 



R' R' R" 

n ab n ad n ac 



1 + 4/ 2 S&' 



B 'ab B 'cd B ac 



J 2 



J (1 + 4/ 2 B'^y 



R' R' R' R"' 
B ac B bc B cd B 



(1 + 4/ 2 £?o ) 3 



(C.21) 



where Sym abcd denotes 1/24 times the sum of all 24 per- as: 
mutations of the indices a, b, c, d, and we note B' Q = B'(0), ~ — 1 



afc' u bcl u ca 



Bq = B "(0) and Sg" = B"'(0). Note that all terms contain- [VVi ~ 2T 2 2^^^ ~ 317^3 2^ 5 ^ 
ing B'(0) can be eliminated by the redefinition B' new (x) = (D 1) 

B (x) — B'(0). The self-consistent equation ( I3.3H then yields, by a similar 

expansion in numbers of sums: 

APPENDIX D: CUMULANT EXPANSION 5B[v ■ «] , / , . , 

- B'[v 2 ab (x) + 2TI 1 



FOR NON-LOCAL EFFECTIVE ACTION 5(v a (x) v b (x)) 

SB[vv] SB[v 



The cumulant expansion can be generalized to study the ef- + 4 / C 2 . 



fective action for non-uniform configurations. The functional 



■nj 



S(v a (y) ■ v b (y)) S(v a (y) ■ v a (y)) 



U[v ■ v] is a functional of the field v 2 b (x) and can be expanded (D.2) 
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APPENDIX E: ANALYSIS FOR ALL m AND 
SHORT RANGE CASE 



9 = 0, i.e. C = (2-d)/2ford < 2or( = Oford = 2, so that 
T m = 4AdT/e does not flow. This is the case studied here. 



To refine the analysis and study the behavior for all values 
of m, let us rewrite from J6. 10b the condition that B remains 
analytic for all to. It reads: 



T>f{m) for all 

{B")- 1 
f(m) = - 



i 

-ilo 



(407 J 2 )- 



2/1 



2/1 



(E.l) 



(E.2) 



where we have also inserted the value of the inverse function, 
with B" ((B")^ 1 (x)) = x, for the power law models. This 
condition is equivalent to the vanishing of the replicon, i.e. it is 
the line where the RS solution of the GVM becomes unstable 
to RSB. 

One can then plot (see Fig. 1101 the function /(to) for the 
three cases defined in the text, LR 0(7) > 0, marginal 0(7) = 
and SR %) < where 0(7) = d - 2 + The LR and 
marginal cases, which correspond to continuous and one step 
marginal RSB solutions, have been discussed in the text and 
there the FRG gives back exactly a m (0) of the MP solution. 
We defer the detailed study of the SR case to further work, 
and give here a few general remarks. 

First one notices on Fig. ^| that solving the FRG equation 
decreasing the mass from infinity one first has the analytic 
solution which coincides with the RS one. For T > T* = 
max m /(to) it remains analytic down to to = 0. For T < T* 
a cusp arises when the left branch of the line is reached. Thus 
despite the reentrance of the analytic solution at small to, 
freezing of the FRG solution has already occurred and it is 
clearly important to understand how to extend the FRG in the 
shaded region. On the other hand, a 1-step solution of the MP 
saddle point equations exists, obtained by varying the free en- 
ergy w.r.t. u c . Its precise boundary depends on the model, but 
it is generally contained within the shaded region (for details 
see illlHllllTS])- An intriguing property of the GVM is 
the simultaneous existence, within the rightmost portion of the 
shaded region near the axis to = 0, of two locally stable so- 
lutions, one RS and the other 1-step RSB. Thus, although the 
line in Fig.^|is the locus of a continuous transition from RS 
to RSB, in this rightmost portion of the to — T diagram, the 1- 
step non-marginal RSB solution appears discontinuously, be- 
fore the line is reached as T is lowered. Work is in progress to 
make contact between FRG and RSB in this SR case, and in 
particular to understand whether there are also two branches 
of solutions of the FRG equation in that region. 



APPENDIX F: FINITE-TEMPERATURE 
FIXED POINTS 

From the RG point of view it is interesting to search for finite 
temperature fixed points (FP) of the FRG equation, especially 
in view of future extensions to finite N (since we know at least 
in some cases, that these persist at finite N). It is convenient 
to use Eq. (17.81 . These FPs exist only in the marginal case 



1. d < 2 

Following the same steps as in Section fVUBI the general so- 



lution of the fixed-point condition md m x(y) = in J7.8I for 
a fixed value of £ > 0, imposing x(yo) = is 



x*{y) = - 

e 



yo 

2C 



e-2C (yo 

■ V v 



2C 
i-2C 



2C(e*o - 1) 



2C 

s-2C 



(Fl) 



Taking a derivative at y = yo, we obtain a self-consistency 
condition for x' Q . One solution is x' — 0, the "zero- 
temperature" fixed point discussed in Section lVlI Bl The other 
one is 



yo(e-2C) 



- 1 



(F.2) 



with the condition that it must be positive. Reinserting in ( IF. Q 
we obtain the final form for the finite-temperature fixed point: 



*, s y yo 

x {y) = i.-x 



1 

y - yo 



(-K yo 
> yo — 

2Ce V y 
yo(e-2C) 



2C 

!-2< 



2C 

!-2< 



T 
2C 



2< 
!-2C 



- -1 



(F.3) 



The term in the second line of dF.3l was not present in the 
"zero-temperature" fixed-point solution ( I7.15> . Note that it 
does not vanish at T = but at the higher temperature T = T c 
such that T m = y °( £ ~ 2< >> , At this point, also x' vanishes and 
the solution b'(x) = —y(x) becomes non-analytic. The fixed- 
point analysis alone does neither fix the value of y , nor T c . 

However we can now explicitly check that this fixed-point 
solution identifies with the analytic solution 17.271 . ( I7.28> 
when setting to — > 0, using T m = T m /(2 — d) and ( = 
I/7 = (2 - d) /2. This identification works only for T > T c , 
and since yo is now fixed by J7.28l >. we can compute T c and 
find that it is given by ( 17.391 1. Below T c the solution freezes at 
to = to c at the zero temperature fixed point. 

2. d= 2 

Let us now solve the fixed point condition md m x(y) = in 
d7.8i for C = 0, imposing x(yo) = 0. One finds: 

y~y° fvo , T m x 'o \ , , , s mA \ 
x {y) = — + — — - My/yo) ■ (R4) 

e V e ex n - 1 / 



Determining x' again one finds, in addition to the solution 

T 

1 -*- n 



x' =0 



-CXn 



1 



yo 



(F.5) 
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— > 
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kf(m) 



9( 7 ) > 




kf(m) \ 








m c\ 


3» 

m 



FIG. 10: The function f(m) defined in Eq. <E.2> . 



0(7) < 



Reinserting one finds finally the finite-temperature fixed 
points 



v - 2/0 



T 

— My/yo) 

67T 



(F.6) 



with e = 2 and T m = T/ 7r. There is thus a line of fixed points 
with ( — in d = 2, parameterized by temperature, yo being 
again undetermined. 

To compare with the solution of the flow equation, we ob- 
viously need to consider a broader class of SR models with 
B'(z) = — gexp(— z/a 2 ). The solution is then: 

x = a 2 ln(y /y) + e- 1 (y-y ) (F.7) 
y = ,gm T /(™ 2 )- £ A- T /(™ 2 > (F.8) 

with e = 2. For small disorder g, and T > T c = eira 2 = 
2na 2 , yo(m) flows to zero as to — > and the solution remains 
analytic. For T < T c the solution develops a cusp when yo 
reaches y — yo(m c ) — ea 2 , i.e. at the Larkin mass: 



(<?/a 2 e) 2 < 1 - T / T c) A t c-t 



(F9) 



Thus only for T = T c the solution reaches for m — * an 
analytic finite-T fixed point associated with £ = 0, of the form 
( IF. 61 . Thus in d = 2 the line of finite-temperature fixed points 
with C — corresponds to the line of critical fixed points as 
the parameter a is varied. 

APPENDIX G: RG FORMULATION OF THE 
RSB SOLUTION 

In this appendix, we derive simple RG equations for the MP 
solution in the full RSB case. This gives a more direct deriva- 
tion of the key-equations and J8.52I . We start from (see 
l lOt and dp . equivalent to (5.4) of MP): 



(G.l) 



a(u) = --B> {2T J (G(k)-G(k,u)) 
Taking a derivative with respect to u yields 



1 = 45" [2T J G(k) - G(k, u) 



k (k 2 + M 



[a](u)) 2 
(G.2) 



for all u such that <r'{u) ^ 0. Using that due to (see (18.5 
equivalent to (5.2) of MP): 



G(k)-G(k,u) = 

we have 

d u {G(k) - G(k,u)) 



u(k 2 + m 2 + [a](u)) 
nl dv 1 



v 2 k 2 



m 2 + [<j](v) 
a'{u) 



(k 2 + m 2 + [<j](u)) 2 



(G.3) 



(G.4) 



On the other hand, we can take a derivative of iG.W with re- 
spect to m 2 : 



d m 2<j{u) = -AB" [ 2T J (G(k) - G(k,uj) 



x / d m 2(G(k)-G(k,u)) 

Jk 

Eliminating B"(. . . ) one finds: 

d m 2a(u) 



(G.5) 



1 



(G.6) 



k (k 2 + TO 2 + [cr](u)) 2 



= - / d m2 (G(k)-G(k,u)). 

Jk 



Taking another derivative with respect to u after using ( IG.4> 
gives 



du 



d m 2<r(u) 



k (k 2 + to 2 + [<j](u)) 2 
lat the derivative; 

ua' (u)d m 2a(u) 



a'(u) 



k (k 2 + TO 2 + [cr](u)) 2 ' 
(G.7) 

Noting that the derivatives of the numerator cancel, we get 



a'(u)(l + d m 2[a](u)) 



k {k 2 + TO 2 + [cr](u)) 3 J k (k 2 + TO 2 + [cr](u)) 3 

(G.8) 

Since for all u, J k (k 2 + m 2 + [a] (u)) 3 ^ and by assump- 
tion a' (u) 7^ 0, RSB reveals its universality in the simple re- 
lation: 



ud m 2<r(u) = 1 + d m 2[a](u) 



(G.9) 
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which, upon another derivation, yields the two "RG equa- 
tions": 



dm 2 



+ [*}(u)) =0 
a 

cr(u) = . 



dm 2 



(G.10) 
(G.ll) 



APPENDIX H: CONVERGENCE TO THE 
FIXED POINT 

Since we have found the solution of the FRG equation for ar- 
bitrary disorder correlations, it is instructive to study the con- 
vergence to the FRG fixed point in the case where the inital 
disorder is not of the simple form J2. 10l > on an explicit exam- 
ple. We start from a disorder correlator which is the superpo- 
sition of two power laws: 



-B'(z) 



9 



C- 



fJ 



(a 2 + z)t (a 2 + z) a 



(H.1) 



with a > 7, s.t. for large z the first term dominates and fixes 
the exponent ( = 2 (i+ 7 ) ■ We w ^ determine the inverse func- 
tion <f>(x) only for C small. One finds 



with C ~ C, and we have defined 

1 1 — a + 7 



7 



(H.2) 



(H.3) 



We have chosen g = (7/e) 7 to simplify all prefactors. 
Inserting into ( 18.421 1. we obtain: 



(HA) 



By multiplying with m e and y Q s , we obtain the equivalent 
formula 



Vo' 



(H.5) 



The left-hand side is plotted on figure^J note that always the 
first exponent is negative, and the second and third are posi- 
tive. Thus the solution for m — or C = is simply yo = 1. 
For non-vanishing C and m, the solution can be obtained 

s — , 

graphically as the intersection of F(j/o) with -CmW^ . 
Note that there is a solution for any C. 

For m — * 0, the approach to y = 1 is obtained by lineariz- 
ing F(y ), with the result 



yo - 1 



jC_ 

1 + 7' 



Ul (1+^)5 f 



(H.6) 



APPENDIX I: PURE O(N) MODELS, 
NON-ANALYTIC EFFECTIVE ACTION 

In this Section we recall the corresponding result for the ef- 
fective action of the generic pure O(N) model at large N. 
One mechanism by which the effective action may becomes 
non-analytic is given on standard example of 4> 4 . 

1 . Self-consistent equation 

The generic O(N) model in dimension d is defined by the 
action: 



s = i 

T 



1 , , „2 1 9 , ,0 fu(x) 2 \ 



Here m is used as a parameter, the bare mass being m 2 



(1.1) 

2 _ 

b — 



2V' (0). For a uniform mode one has T[u] = L T[v] in 



(1.2) 



terms of the rescaled field v — u/yN. One defines 

f[v] = imV + V(v 2 ) = W(v 2 ) . 
Similarly, one defines 



W(z) = -m 2 z + V(z) 



(1-3) 



whether absorbing or not the mass into the (bare or renormal- 
ized) potential. Again, for m = 00, one has V = V. The 
same method as in Section dlllt yields the saddle point equa- 
tions for infinite N: 




V'(v 2 ) = V'(v 2 + G(v)) 



G(v) 



T 



q q 2 + m 2 + 2V'(v 2 +G(v)) 



(1.4) 
(1.5) 



More details, a graphical derivation, and the 1/N expansion 
are given in [47 ] . A condition for the stability of the theory is 
that: 



2W'(v 2 ) := m 2 + 2V'(v 2 ) > for all v 2 

2. Solution and FRG equation 

Let us start from the form 



(1.6) 



FIG. 11: The function <H.5> . describing the approach to the fixed 
point in presence of an additional bare power law tail. 



W'(x) = W [x + T 



1 



q q 2 + 2W'(x) 



(1.7) 
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The self -consistent solution of this equation is formally: 



x=(W'r 1 {y)-T 
V = W'(x) . 



1 



, q 2 + 2y 



(1.8) 
(1-9) 



and y > 0. 

Let us write the associated FRG equation. One has 



md m W'(x) = -rn 2 - 2TW" [x + T 



1 



x ^— md m W'(x) 



W"(x) = W" [x + T 



q q 2 + 2W'(x) J 
1 

q (q 2 + 2W'{x)) 2 

(1. 10) 

1 \ 



q q 2 + 2W'(x) 



x 1 - 2TW"(x) 



, (q 2 + 2W'(x)) 2 



(I.H) 



Thus: 



-md m W'{x) = 2Tm 2 W"(x) 



(q 2 + 2W'{x)f 



(1.12) 



For d < 4 taking A to infinity, this becomes 



-md m W'(x) = 2T—m 2 W"{x)(2W'(x))-^ 2 - m 2 



(1.13) 



3. </> 4 -theory and non-analytic behavior 

For the 4 theory V(x) = | (x — l) 2 this reads: 

W'(x) = \m 2 + |(z - 1) + ^ / A ^ 

K ' 2 2 V 2 J q q 2 + 2W'(x) 

, x 2 m 2 r A 1 

^(2/) = -y + 1 r / 

3 9 Jq 



q 2 + 2y 



(1.14) 
(1.15) 



for y > yo with x(yo) = 0. At T = 0, as m is decreased, 
there is a transition in any <i when m 2 = m 2 — g vanishes. 
For d > 2, the transition persists for T < T c , and occurs 
when m 2 — g + gT f ^ vanishes, with the standard result: 



(1.16) 



Tel ? = 1, 



which depends strongly on the UV-cutoff A. yo vanishes at the 
transition, and in the ordered phase the effective action has a 
non-analytic form. In addition of the branch ( II. 1 51 for y > 0, 

x > x c = 1 — — g T the function x(y) has a branch 

y = for < x < x c , where x c is the order parameter. 

Exactly at T = T c we should recover that the effective ac- 
tion exhibits the standard power-law non-analyticity T[M] = 




FIG. 12: Relation between x and y in scalar field theory: above, at, 
and below the critical temperature. 



M 1+s . Indeed, from the self-consistent solution dI.15> . sub- 
tracting the same with x = x c and y = one gets 



9 r A / 2 y i i 

x(y) = -y+(2yr- 2 y 2 T c 1 



fJ 



k 2 k 2 + 1 



(1.17) 



for 2 < d < 4. This corresponds to S = (2 + d)/{2 - d). It 
can be recovered by solving the FRG equation. One can look 
for fixed-point solutions of the form 



W'(x) = m a f(m' j x) . 



(1.18) 



If one wants the two first terms to dominate and to scale in the 
same way, one needs a < 2 and f3 = ae/2 — 2. For all three 
terms to scale the same way one needs a — 2, j3 = 2 — d. 
Inserting dl. 1 8t into dl. 1 3t yields 

1 - af(z) - Pzf'(z) = 2T^f'(z)(2f(z)r^ 2 . (1.19) 

Again this can be transformed into a linear RG equation for 

*(/): 



(1 - af)z'(f) - pz{}) = 2T^(2/)-/ 2 



(1.20) 



The solutions of the above equation with a = 2, (3 = 2 — d 
are: 

= _2T^(2/ - l)( d - 2 )/ 2 f t^' 2 (2t - iy d / 2 



A particular solution is 

z(/) = (2/-l)( rf - 2 )/ 2 
W'{x) = ?^-+ x 2 /( d - 2 ) . 



(1.21) 

(1.22) 
(1.23) 



In the limit of zero mass this yields V'{x) = x 2 ^ d ^ 2 \ 

One can also pursue the RG approach in the ordered phase, 
as is done usually in the form of a non-linear sigma model, 
and deal with a non-analytic effective action. 
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Although the mechanism for the disordered systems stud- 
ied in the main text seems to be different from 4 -models, it 
raises the question of the meaning of the non-analyticity in 
the disordered problem. Is it the signature that we are dealing 
with a glass phase, where a symmetry has been broken? We 



know that for infinite N, this is also accompanied by RSB, 
but this does not have to be so in general, i.e. the cusp can 
arise without RSB, just from localization (single ground state 
dominance). 
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